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Euler Equation on a Rotating Surface 

Michael Taylor * 


With an Appendix by Jeremy Marzuola and Michael Taylor 

Abstract 

We study 2D Euler equations on a rotating surface, subject to the 
effect of the Coriolis force, with an emphasis on surfaces of revolution. 

We bring in conservation laws that yield long time estimates on solu¬ 
tions to the Euler equation, and examine ways in which the solutions 
behave like zonal fields, building on work of B. Cheng and A. Mahalov, 
examining how such 2D Euler equations can account for the observed 
band structure of rapidly rotating planets. Specific results include 
both an analysis of time averages of solutions and a study of stability 
of stationary zonal fields. The latter study includes both analytical 
and numerical work. 
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1 Introduction 


Let M = do be a surface in rotating about the xs-axis at constant an¬ 
gular velocity ui = —0/2. A natural class of such bodies would be those that 
are rotationally symmetric about the xs-axis, and we will eventually settle 
into the study of that class, but initially we will not make that assumption. 
We will assume M is diffeomorphic to the standard unit sphere S'^. We aim 
to study 2D incompressible Euler flows on M. 

The approach of Rossby to the effect of the Coriolis force on flows on 
M, described on p. 21 of yields the Euler equation 

dxi 

—+ VuU = Qx{x)Ju — Vp, divu = 0, (1.0.1) 

where 

X(x) = 63 • u{x), ( 1 . 0 . 2 ) 

zz(x) being the unit outward pointing normal to M at x. Here u is the flow 
velocity, a tangent vector field to M, and J : T^M —)• T^M is counterclock¬ 
wise rotation by 90°. In case M = S'^, we have x(x) = X 3 . For more general 
M that are rotationally symmetric about the X 3 -axis and that have positive 
Gauss curvature, we have x(^) = 

Bringing in the 1-form u, arising from u via the isomorphism TxM ~ 
T*M determined by the metric tensor on M, we can rewrite (11.0. 1|) as 


--—h = Dx * u — dp, (5u = 0. (1.0.3) 

at 

We can eliminate p from (I1.0.3P via the Leray projection P, the orthogonal 
projection of L^(M, A^) onto the subspace where 5u = 0. We get 


du 

— -h PVuU = QBu, u = Pu, 

dt 


(1.0.4) 
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where 


Bu = P{x*Pu). (1.0.5) 

We mention a few essential properties of B, which will facilitate the 
analysis of (11.0.41) . First, one can deduce from the Hodge decomposition 
that, on 1-forms, 

P = -5A2^d, (1.0.6) 

where denotes the inverse of the Hodge Laplacian on 2-forms, dehned 
to annihilate the area form. Thus 


Bu = —5^2 ^d{x * Pu) 

= —6^2^{dx A *Pu)^ 


(1.0.7) 


since d * Pu = 0. We deduce that 


H is a compact, skew-adjoint operator on L^(M, A^). (1.0.8) 

In fact, B G OPS~^{M), i.e., H is a pseudodifferential operator of order 
— 1. The skew adjointness is a direct consequence of the formula (jl.O.Sp . the 
skew adjointness of the Hodge star operator *, and the commutativity of * 
and multiplication by x- Further results on B can be found in ^1.11 

Our interest in the equation (ll.O.ip was stimulated by the recent paper 
[6] of B. Cheng and A. Mahalov, investigating the case where M is the 
standard sphere 5^. That paper took (jl.O.ip as a model of the behavior of 
the atmosphere of a rotating planet, and investigated how it might account 
for observed band structure, particularly on rapidly rotating planets, such 
as Jupiter. This involved a study of zonal flows, i.e., velocity helds of the 
form JV/, where / = /{x^) is a zonal function. The paper looks at time 
averages, 

1 

As,tu = f I u{t) dt, (1.0.9) 

J s 

for a solution u to (jl.O.ip . The main result (Theorem 1.1 of [6]) is that, 
if uq G k > 3, divuo = 0, there exists Tq > 0, independent of 

H, such that (11.0.Ih has a unique solution for t G [0, To/||uo||j|/fc], satisfying 
u(0) = Uq, and, for0<S'<S'-|-T< To/||uo||//fe, one has 

11(1 - U)As,tu\\h^-^ = 0(|lir^), (1.0.10) 

where H is a projection of the space of divergence-free velocity helds on 
onto the space of zonal helds. 
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In the present paper, we push the study of (|1.0.ip in the following direc¬ 
tions. 

(A) Investigate a larger class of rotating bodies, beyond the class of rotating 
spheres. 

(B) Establish estimates on As,tu that are uniform in S', T S (0, oo), without 
restrictions on their size. 

(C) Investigate another way that large |n| enhances band formation, namely 
by enhancing the stability of zonal fields as stationary solutions to (|1.0.ip . 

These are natural directions to pursue. Rapidly rotating planets are 
flattened at the poles and bulging at the equator. Also, a planet like Jupiter 
has been rotating for a very long time. Of course, of major interest to us is 
the set of interesting new mathematical challenges that arise in addressing 
these issues. 

We proceed as follows. In ^we produce basic existence results, starting 
with short time existence in Results of ^ apply to any surface M 

diffeomorphic to 5^, with no symmetry hypothesis on the geometry. To go 
from short time to long time existence, we derive in l)2.2l an equation for the 
vorticity w = rot u, namely 

^{w-^x)+ Vu{w-^x) = 0. ( 1 . 0 . 11 ) 

This is a conservation law, yielding a uniform bound on ||r(;(t)||ioo on any 
time interval on which (jl.O.ip has a sufficiently smooth solution. Using this, 
we adapt the classical Beale-Kato-Majda argument [2] to establish existence 
for all t of a solution to (ll.O.ip . provided tt(0) = uq is divergence-free and 
belongs to H^{M) for some s > 2. This is accompanied by the estimate 

\\u{t)\\Hs <C\\u{0)\\hs exp exp(^C's(||r(;(0)||Loo -hCiriDty (1.0.12) 

In P we specialize to the class of smooth compact surfaces M C that 
are invariant under the group of rotations about the xs-axis, and that in 
addition have positive Gauss curvature everywhere. This hypothesis will be 
in effect for the rest of the paper. As already mentioned, this symmetry 
hypothesis implies x = xi^s) (U-O-ID . In 113.11 we show that when / is 
a zonal function, the associated zonal vector field u = f is a stationary 
solution to (jl.O.ip . We also give examples of stationary solutions that are 
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not zonal fields. In ^3.21 we study time averages of the form (11.0.91) and 
establish estimates of the form 

11(7 - Il)As^Tu\\H-3,qm 

- ^ ^ C||n(0)f-^(||^/;(0)|U^ + 2|fl|)''}, 

(1.0.13) 

for 0 < 1 < 0, with q{9) = 1/(1 — 9). This is valid for all 5, T G (0,oo). It 
should be expected that the norm on the left side of (|1.0.13p is weaker than 
that in (jl.O.lOp . In any case, having a weak norm seems consistent with the 
appearance of complicated eddies within the bands of a planet like Jupiter. 
We proceed in ll3.3l to derive an additional conservation law, of the form 
^ /» 

— J ^{x)w{t,x)dS{x) = 0, (1.0.14) 

M 

for solutions to (II.0.11) on our radially symmetric domain. We discuss com¬ 
putations of y and ^ in 1)3.41 and technical smoothness results in 1)3.51 

In ^ we take up the study of stability of stationary zonal solutions to 
Ol.O.ip . assuming M is radially symmetric and has positive Gauss curvature. 
In ^4.11 we apply an Arnold-type approach, and deduce that a sufficient 
condition for stability in H^{M) is that 

u;(^) — nx(0 is strictly monotone in (1.0.15) 

where w = Af is the vorticity. In ^4.2l we study the linearization at a steady 
zonal solution of (jl.O.ip , or more precisely of the vorticity equation (jl.O.lip , 
obtaining a linear equation of the form dC,/dt = r(/. We establish a version 
of the Rayleigh criterion, namely, if the spectrum of T is not contained in 
the imaginary axis, then 

w'{^) — ^x'iO must change sign. (1.0.16) 

Note that (11.0.151) and (11.0.161) are almost perfectly complementary. Never¬ 
theless, the criterion 01.0.16)) . while necessary for failure of SpecT C iM, is 
not sufficient. This matter is discussed in ^ 

In ^we specialize to M = 5^ and perform some specific computations, 
taking f{x) = cP^{x 3 ), v = 2,3,4. We make use of classical results on 
spherical harmonics to produce infinite matrix representations of the linear 
operator T. We present some analytical results for v = 2 and some numerical 
results for i/ = 3 and 4, indicating how the Rayleigh-type criterion (jl.O.lOp is 
not definitive as a criterion for linear instability. We also discuss the extent 
to which stability seems to depend monotonically on 71 (or, sometimes, not). 
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1.1 Further properties of the operator B 


The operator i?, defined in (ll.U.5|) . arose in the form (I1.U.4I) of the Euler 
equation. By (jl.O.ip . we have 

B £0PS-^{M), (1.1.1) 

the class of pseudodifferential operators on order —1 on M. We record some 
other properties of B, which will be useful later on. 

Since 6u = 0 on M ^ u = *df for a scalar function / (known as the 
stream function), uniquely determined up to an additive constant, it is useful 
to compute 

B{*df) = dA^Hdx^df). ( 1 . 1 . 2 ) 

We have (with a denoting the area form on M) 

dx^df = - * *dx A df 
= df A *{*dx) 

= {df,*dx)a (1.1.3) 

= {df,JVx)a 

= *Zf, 

with the vector field Z given by 



^ = JVx- 

(1.1.4) 

Note that 

The formula (|1.1.2p vields 

div Z = 0. 

(1.1.5) 


B{*df) = * Zf 

= *dAQ^Zf. 

(1.1.6) 


Note that (jl.l.Sp implies that Z is skew-adjoint and that ZfdS = 0. We 
see from (jl.1.61) that 

V n KerB = {*df : f G H\M), Zf = 0}, (1.1.7) 


where 

V = {uG L^{M, A^) ; 6u = 0}. (1.1.8) 

When M = 5^, we have the following result, observed in [6]. 
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Proposition 1.1.1 If M = S'^, then B commutes with Ai, the Hodge 
Laplacian on 1-forms. 


Proof. In such a case, we have (1.1.4) with x(x) = X3, hence Z = A3, the 
vector field generating the 27r-periodic rotation about the xs-axis. Since the 
flow generated by A3 consists of isometries on S'^, A3 and Aq commute. 
Then (by (ll.l.Bp i 

AiB{*df) = *dAoAo ^Ag/ 


= ^dA^^XsAof 

= B{*dAof) 

= BAi * df. 


(1.1.9) 


□ 


2 Basic existence results 

Here we establish existence of solutions 12(0) to (|1.0.3p . given uq G 
s > 2, and 6uo = 0, and we produce estimates on such solutions. We 
begin in ^2.11 with short time existence results. In preparation for long time 
existence results, we derive a vorticity equation in ^2.21 We show that if u 
solves (|1.0.ip and w{t) = rotii(t), then 

^(w - Hx) + Vn(u; - Hx) = 0. (2.0.1) 

This is a conservation law. We use it, together with a method pioneered by 
[2] , in ^2.31 to establish long time existence of solutions to (I1.0.3P . We show 
these solutions satisfy the estimate 

\\u{t)\\H^ < C||12(0)|||^s exp exp (^C^(||r(;(0)||Loc + C|H|)t^. (2.0.2) 

The path from (j2.0.ip to (I2.0.2P passes through the estimate 

||M(t)||(,i.oc < C\\w{t)\\L^ < C(||u;(0)||loo +2|H|), (2.0.3) 

which will see futher use in ^ Here, 

= {u£ C{M) : Vu G bmo(M)}. (2.0.4) 
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2.1 Short time existence 


Our approach to the short time solvability of 111. U. II) , or equivalently (ll.U.ip , 
i.e., 


du 



— + PVuU = QBu, 
ot 

u = Pu, 

(2.1.1) 

with initial data 



u(0) = ho G 

o' 

II 

O 

(2.1.2) 

is to take a mollifier = (/j(eAi), real valued and in 

C'^(M), with (/?(0) = 


1 (Ai the Hodge Laplacian on 1-forms), and solve 
du 

PJeVu.JeUe = ^JeSJ^Ue, 
Ot 


Compare the treatment in §2, Chapter 17, of [12] (for fl = 0). Given e > 0, 
the short-time solvability of (|2.1.3p is elementary, since (12.1.31) is essentially 
a finite system of ODEs. The goal is to obtain estimates of Ue{t) in 
for t in some interval, independent of e, and pass to the limit. 

To start, we have 


2 — (T’JeVug JjUg, Ug)-|-H(JjHJg'Uj, Ug). (2.1.4) 

As in m, the first term on the right is 0 (cf. (2.3)-(2.5) in Chapter 17 of 
|12j). By (ll.O.Sp . so is the second term on the right side of (I2.1.4|) . Hence 


Us{t)\\L^ = ||J£Uo||l2. 


(2.1.5) 


This is enough to guarantee global existence of solutions to ()2.1.3p . for each 

e > 0. 

To estimate higher-order Sobolev norms, we bring in 

A = (-Ai)^/2^ ||h||H= = ||A^u||i2. (2.1.6) 


Then 

= -{A^PJe^uJeUe^A^Ue) 

+ n{A^J,BJ,u„A^u,). 

Now, by (ll.O.Sp . 


{A^JeBJ,U„A^U,) = {BJ,A^U,,J,A^Ue) 

+ ([A^H]J,u„J,A%) 


(2.1.7) 


( 2 . 1 . 8 ) 
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Furthermore, since has scalar principal symbol, 

[A^,B]eOPS^-^{M). (2.1.9) 

It follows that the second term on the right side of (I2.1.7P is 

<C\n\-\\Ue\\jjs-i. (2.1.10) 

We can write the first term on the right side of ()2.1.7p as 

- {A^PJeVueJeUe^A^Ue) = -{A^VuJeUe,A^JeUe). ( 2 . 1 . 11 ) 

In order to make use of the identity {Vu^v,v) = 0, we need to analyze the 
commutator [A^,Vu^]. We claim that 

||[^*, V„Jj££(e||i2 < C\\Ue{t)\\ci\\u^{t)\\H-, (2.1.12) 

with C independent of e. If s = 2k is a positive even integer, this is a Moser 
estimate, as in (2.11)-(2.13) of [12], Chapter 17. For general real s > 0, this 
is a Kato-Ponce estimate, established in [7] in the Euclidean space setting, 
and in greater generality (directly applicable to the setting here) in §3.6 of 

m- 

In more detail, the KP-estimate gives, for s > 0, 

— fA^v\\L2 < C'||/||(7i||u||jys-i + C'||/||//s||u||loo. (2.1.13) 

We take v = Xu, where X is a hrst-order differential operator, and write 

A^fXu) - fX{A^u) = A\fXu) - fA^Xu) + f[A^,X]u. (2.1.14) 

Then (|2.1.13p applies to estimate the first two terms on the right side of 
(|2.1.14p . and the L^-norm of the last term is bounded by CII/Hloo|| ri||. 
Thus 


||[7l^/X]n||i2 <C'||/||ci||n||H^+C||/||^,.||n||ci, (2.1.15) 

which in turn yields (12.1.121) . 

From (|2.1.121) , we bound the absolute value of (12.1.lip hy C\\u^{t)\\ci\\u^{t)\\'jjs- 
Together with (I2.1.10p . this gives 

< C||i,(()||c.||S.(t)ll?,. +C|f!| . (2.1.16) 
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On the 2D manifold M, ||u||( 7 i < C's||'u||h», as long as s > 2, so (12.1.161) 
implies 


-^\\Ue{'t)\\‘H3 < C\\Ue{'t)\\]j3 +C|D| • \\Ue{t)\\‘jj3- 
By Gronwall’s inequality, we have, for t > 0, 






< y(t), 


where y{t) solves 


^ + |D|y), y(0) = ||u£(0)||^3. 


(2.1.17) 

(2.1.18) 

(2.1.19) 


In particular, {u^{t) : 0 < t < T^} is uniformly bounded, in inde¬ 

pendent of e S (0,1], as long as 


Tm < 


hf 

^ Jvi 


dy 


( 2 . 1 . 20 ) 


C Jy{o) + my' 

For a more explicit (though cruder) upper bound, we can say that 

||il.(i)|li3 + l<2(t), (2.1.21) 

where z{t) solves 


dy 

— = (7(1 -I- |D|)2;^/^, 2;(0) = ||u£(0)|||^3 + 1- 


( 2 . 1 . 22 ) 


Explicit integration gives 


z(t) = ^z(O) — (7i(l-|-|D|)t^ , for 0 < t < (7]^ ^2;(0) ^/^(1-|-|D|) 

(2.1.23) 

Consequently, 

||tte:(t)||ci < Nn{t) = ^ _ forO<t<Tm, (2.1.24) 

with 

Ci(l + |D|)’ ^ = (ll^o||^3 + l)-^/". (2.1.25) 

This plugs into (j2.1.16|) to yield 


Ue{t)\\Hs < CNn{t)\\Ue{t)\\Hs -h (7|D| • \\Us{t)\\jjs-i 

<C{Nn{t) + \mMt)fHs, 


(2.1.26) 


10 









for t € [0, Tm), which in turn yields 

\\us{t)\\Hs < \\us{0)\\hs eiQp C f {Nn{s)+ \n\)ds, for t e [0,Tm), 

Jo 

(2.1.27) 

an estimate that is uniform in e G (0,1]. 

With these uniform estimates in hand, one can apply standard tech¬ 
niques, discussed in Chapter 17 of [12], to obtain a solution u{t) to (I1.U.3I) 
in C{[0,Tm), (M)), given initial data uq G H^{M) (satisfying 6uo = 0) 

as long as s > 3. Here Tm is as in (|2.1.25|1 . and the solution u{t) satisfies 
an estimate parallel to (I2.1.27p . Also, estimates parallel to those produced 
above establish uniqueness of the solution u{t) and continuous dependence 
on the initial data uq. 

Remark. One could replace in (I2.1.17P by for any sq > 2, and have 
a local existence result for initial data uq G H^{M) for any s > sq. 

Improved estimates for M = S^ 

The estimates (I2.1.25P and (|2.1.271) for the existence time and size of 
solutions to (I1.0.3P exhibit a dependence on |R|. It was observed in |6| that 
one has estimates independent of R when M = S'^ is the standard sphere. 
We note the changes in the arguments above that yield this. 

The key modification arises in the estimate (I2.1.10p on the second term 
on the right side of (I2.1.7p . If M = S'^, then B commutes with Ai (cf. Propo¬ 
sition [LTT]) , hence with and J^, so 

{A^JeBJ.u,, A^u,) = {BA^J,u,,A^JSe) = 0, (2.1.28) 

the latter identity by the skew adjointness of B. Hence (12.1.101) is replaced 
by 0, and (j2.1.161) is improved to 

j^WMOWh. < C'||«,(()||c.||%(*)||l,., (2.1.29) 

provided M = S'^ and s > 3. In this case, Gronwall’s inequality yields 
(12.1.181) where y{t) solves 

^ = y(0) = ||n,(0)||^3. (2.1.30) 

This has the explicit solution 

2/(f) = (y(0)-i/2 _ CR)-2, for 0 < t < C^^y{0)-^/\ (2.1.31) 
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Consequently, (I2.1.24P is improved 


\\ue{t)\\c^ < N (t) = for 0 < t < Tm, (2.1.32) 

with 

Tm = C{^A, A = \\uo\\-l, (2.1.33) 

and ()2.1.27p is improved to 

N{s)ds, for ts[0,Tm), (2.1.34) 

given s > 3, an estimate that is uniform in both £ G (0,1] and 11 G M. From 
here, arguments parallel to those indicated above give a unique solution to 
(jl.0.3p . with initial data ho G for s > 3, on t G [0,Tm), satisfying 

an estimate parallel to (I2.1.34p . This result is similar to Theorem 5.1 in 
[6], except that here (thanks to Moser-type estimates) the t interval is in¬ 
dependent of s > 3 (and s is not required to be an integer, and also we can 
actually fix sq > 2 and take s > sq). 

2.2 Vorticity equation 

The Euler equation (I1.0.3P can be rewritten as 

du /I \ 

— + CuU = u + d^-\u\^ - pj , 5u = 0, (2-2.1) 

where C is the Lie derivative. This follows from the general identity 

V„h = —-d|u|^. (2.2.2) 

Compare [l2]. Chapter 17, §1. We obtain an equation for the vorticity rc, 
given by 

du = w = wa, w = r:otu, (2.2.3) 

where a is the area form on M, by applying the exterior derivative to (j2.2.ip : 

dw 

— -h LuW = nd{x * u) 

= Q{dxA*u) 

= ^{dx,u)a, 


it)\\H. < II«£(0)||h» exp C 


i‘ 
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hence (since CuOl = 0), we have the vorticity equation 


dw 

— + VuW 

ot 


^{dx,u) 

^VuX- 


(2.2.5) 


We can rewrite (j2.2.5j) as 

^{w-^x)+ Vu{w-^x) = 0 , ( 2 . 2 . 6 ) 

which is a conservation law. 

It is useful to know that we can reverse the path from (j2.2.ip to (j2.2.4p . 

Proposition 2.2.1 Assume 5u = 0 and set w = du. If w satisfies (j2.2.4jl . 
then u satisfies (I2.2.ip . 

Proof. For such u, the Hodge decomposition on M allows us to write 

du ~ 

— +CuU-nx*u = dF + G, (2.2.7) 

ot 

where G is a 1-form on M (for each t) satisfying 

5G = 0. (2.2.8) 

Applying the exterior derivative to (I2.2.7p yields 

d " 

+ C-uW = H(Vttx)® + dG. (2.2.9) 

If (j2.2.4p holds, we deduce that 

dG = 0, (2.2.10) 

which, in concert with ()2.2.8p . implies G = 0, since the hypothesis that M 
is diffeomorphic to 5^ implies 77^(M, M) = 0. □ 

Also the identity = 0 allows us to write 

6u = 0 => u = *df, hence u = JVf, (2.2.11) 


with scalar / (the stream function) determined uniquely up to an additive 
constant, which we can specify uniquely by requiring 


















Note that 

w = A/, (2.2.13) 

and we have 

u = *d/S.~^w, (2.2.14) 

with A~^ dehned on scalar functions to annihilate constants and have range 
satisfying ()2.2.12l) . We can rewrite the vorticity equation ()2.2.5p as 


dw 

'm 


+ {jvf,v{w-nx)) = o. 


(2.2.15) 


2.3 Long time existence 

As seen in 1)2.11 if uq G s > 3 (or even s > 2) and Suq = 0, then 

(|1.0.3p has a solution u G C{[0,Tq), satisfying 

< C|l«(()llc. IlilWIlH. + Clfil ■ ||u(()l&.-. „ 3 

< C\\u{t)\\ci\\u{t)\\H‘ + • \\'d{t)\\Hs, 

for some Tq > 0. The analysis behind short time existence shows that if 
[0, Tq) C M"*" is the maximal interval of existence for u, with such regularity, 
and To < oo, then ||ti(t)||//s cannot remain bounded as t Tq. 

Our goal here is to demonstrate global existence of such a solution. We 
use the method of [2] to obtain such long time existence. (An alternative 
approach could proceed as in the analysis in m-) A key ingredient is the 
vorticity equation (|2.2.6p . which is a conservation law. It implies that, for 
all t G [0, Tq), 

||t(;(t) - Ox||l°o = ||t(;(0) - 0,x\\l^, (2.3.2) 

where w{t) = rot u{t). It follows that 



?i^(i) < Kl>(0) + 2 11 , 

(2.3.3) 

since, by 

(ll.0.2h. Y r°° = 1- Now u+fl foo does not bound 

ll^tWIIcT but. 

since 

u{t) =A:[^6 € 0PS-\M), 

(2.3.4) 

we have 

\\u{t)\\cl < ||u(t)||(,i,oc < C\\w{t)\\L<^, 

(2.3.5) 

where 

= {u£ C{M) : Vu G bmo(M)}, 

(2.3.6) 
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and Cl{M) is a Zygmund space. A variant of the analysis of [2] (See [13] 
Appendix B) gives 


|ft||ci < C\\u\\ci (l + log ), ?’>0. 


Hence 


||ft||ci < C'll«llci(i + iog'^ II''^IIh=)> 

provided s > 2. Plugging into (I2.3.ip . we get, for s > 2, 


^1111(^)11^3 < C's(lk(0)||Loo + C'|fl|)(l + log+ \\u 


Now, with 
(j2.3.9p says 


^ < ^^^(l + log+ y{t))y{t), 


so 


rvit) 

Jv(0) 


drj 


Now, for y > e, 


y{0) r?(l + log+r/) 
y drj 


< CsAt. 


[ 


ylog rj 


= log log y, 


so 


(2.3.7) 

(2.3.8) 


Hs)\\u[t)\rHs. (2.3.9) 
yit) = Mt)\\h, A= ||u;(0)||l-+C| fl|, (2.3.10) 


(2.3.11) 

(2.3.12) 

(2.3.13) 

(2.3.14) 


y{t) < if y{0) = e. 

From this we can deduce that 

\\u{t)\\Hs < C||it(0)|||^a exp exp(^C5(||r(;(0)||Loo +C|n|)ty (2.3.15) 

This estimate implies that 1112(1)is bounded on [0, Tq) for all Tq < oo, so 
we have global existence, with the global estimate p2.3.15p . 


Remark. As seen in ^2.11 when M = S'^ one has an improvement on (|2.3.ip . 
namely, the term on the right side containing |R| can be dropped. However, 
the term containing |R| in p2.3.3p persists, so one does not get a signihcant 
improvement on p2.3.9p . or on (I2.3.15p . in this case. 
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3 Bodies with rotational symmetry 

Here we assume M C is invariant under the group 

( cos s — sin s \ 

sins coss (3.0.1) 

1/ 

of rotations about the xs-axis. We also assume that M is diffeomorphic 
to and has positive Gauss curvature everywhere. We consider special 
properties of the Euler equation (jl.O.ip under this hypothesis. Note that, if 
X is given by (ll.0.2p . then 

Xsy = 0, (3.0.2) 

where X3 denotes the vector field on M generating the flow (13.0.11) . In fact, 
we can write 

X{x) = x(x3), (3.0.3) 

where on the right side x G C'°°([—a, a]), assuming 

X 3 : M ^ [—a, a], a = rnnx X3, —a = nun X3, (3.0.4) 

which can be arranged by a translation. In such a case, the vector field 

Z = JVx, (3.0.5) 

arising in (ll.l.dp . is parallel to X 3 . In fact, 

Z = 4>X3, 4>(x) = -x'ixs)- (3.0.6) 

If M = 5^, then xix) = X 3 , and Z = —X 3 . 

In 1)3.11 we study stationary solutions to (ll.O.ip . i.e., solutions that are 
independent of t. We show that if / G C°°{M) is a zonal function, i.e., 
X^f = 0, then the associated divergence-free vector field u = JX f (which 
we call a zonal field) is a stationary solution to (ll.O.ip . for all H. We also 
give examples of stationary solutions that are not zonal fields. 

In 1)3.21 we return to time-dependent solutions and study time averages 

As,tu = I u(t) dt, (3.0.7) 

'j s 

where u, solving (ll.0.3p . is the 1-form counterpart to the vector field u, 
solving (jl.O.lj) . We construct a projection H from forms satisfying 5u = 0 
onto the subspace of zonal forms, and produce estimates on 

II(/ - II)As,tu\\h-3,<i , (3.0.8) 
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for 1 < g < oo, involving negative powers of |n| (see ()3.2.36p l. Our interest 
in such estimates was stimulated by the paper [6], which produced estimates 
on ||(/ — YV)As,Tu\\Hk ,2 for positive k, valid however over a limited range of 
S, T. Our goal was to produce uniform estimates, valid for all time. One 
key to this was to use the long-time existence and estimates, from 112.31 in 
place of short-time existence and estimates, from ^2.11 Also, the results 
of were derived for a rotating sphere. Since fast rotating planets have 
noticeable bulges at the equator, we were motivated to treat more general 
rotationally symmetric cases M. 

In ^3.31 we produce another conservation law. Namely, with ^ S C°°{M) 
satisfying 

^3 = (3.0.9) 

if u satishes (11.0. ip and w = rot u is the associated vorticity, then 

J ^{x)w{t,x) dS{x) (3.0.10) 

M 

is independent of t. Note that M = ^ ^(x) = xi^) = ^ 3 . Such a 

conservation law appears in [CaM] in the special case M = S'^. We derive 
it here for a similar reason as [CaM], as a tool to use in an Arnold-type 
analysis of stability of stationary solutions to (|1.0.ip : see @ 

In ^3.41 we discuss computations of x and hrst for a general surface of 
revolution 

xl + xl = r{x 3 f, (3.0.11) 

and then, more explicitly, for ellipsoids of revolution 

xi + xl + = 1. (3.0.12) 

Section [331 establishes smoothness of various functions, such as xixs) and 
^(xa), making essential use of the positive Gauss curvature assumption. 

3.1 Stationary solutions 

A stationary solution to (jl.O.ip is one for which dujdt = 0. In such a case, 
w = rotn satisfies (|2.2.5I) with dw/dt = 0. Hence, by ()2.2.15p . 

(JV/,V(u;-Hx)) = 0, (3.1.1) 

where 

w = A/, u = *df, (3.1.2) 
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which determines / uniquely, up to an additive constant. The equation 
(I3.1.ip is equivalent to 

V(A/- II V/. (3.1.3) 

By Proposition 12.2.11 whenever / satisfies (I3.1.3p . which implies (j3.1.ip . 
then u, defined by (I3.1.2p . is a stationary solution to (jl.0.3p . This gives 
the following class of stationary solutions. We say / G C°°(M) is a zonal 
function if X 3 = 0, where the vector field X 3 generates 27r-periodic rotation 
about the xs-axis. Then we say u = JV/ is a zonal velocity field. 

Proposition 3.1.1 Assume M <Z M? is a smooth, compact surface, with 
positive Gauss curvature, and radially symmeric about the x^-axis. If f ^ 
C°°{M) is a zonal function, then u = JV/ is a stationary solution to (jl.O.ljl . 
for all Q. 

Proof. Under our hypothesis, we have x = f = and w = 

w{xs), so (13.1.31) holds. □ 


Corollary 3.1.2 With M as in Proposition 

each u £ V r\ Keic B is a stationary solution to (ll.0.3p . (3.1.4) 

Proof. Recall that V n KerR is given by (|1.1.7I) . with Z = JVx, as in 
(|1.1.4p . The geometrical hypothesis on M made above implies 


Z = 4>X3, (3.1.5) 

for some nowhere vanishing G C°°{M), which yields (I3.1.4p . □ 

While our study of stationary solutions to the Euler equation will focus 
on zonal functions, we mention that there are stationary solutions to (11.0.31) 
that are not of the form (j3.1.4j) . We give examples when M = 5'^, the 
standard sphere. To get started, note that (I3.1.3P holds whenever there is a 
smooth : M —^ M such that 

+ ^3:3 (given M = 5^). (3.1.6) 

We will apply this with V^(/) = —Xkf, where Xk is chosen from 

Spec(-A) = {Xk = k‘^ + k : k = 0,l,2,3,...}. (3.1.7) 
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Note that X 3 is an eigenfunction of —A with eigenvalue Ai = 2. Thus we 
assume k >2. Then (I3.1.6P becomes 


(A + Afc)/ = Qxs- 


(3.1.8) 


As long as A^ / 2, (|3.1.8p has solutions, and the general solution is of the 
form 



Thus g is the restriction to of a harmonic polynomial, homogeneous of 
degree k. For example, we can take 


/c = 2, Afc = 6, gkix) = xl-xl, 
k = 3, Xk = 12, gkix) = Re(xi + ix 2 f 


(3.1.10) 


etc. For such / as in (|3.1.9p . we have 



(3.1.11) 


as a stationary solution to (jl.O.ip . 

3.2 Time averages of solutions 

As before, M C is a smooth compact surface of positive Gauss curvature 
that is radially symmetric about the xa-axis. We take u to be a smooth 
solution to the Euler equation (11.0.lb . so u solves (ll.O.Op . or equivalently 


(ITTO . i.e. 



(3.2.1) 


Given S,T £ (0,00), we want to investigate the time-averaged field 



(3.2.2) 


In particular, we investigate the extent to which it can be shown that As,tu 
is close to a zonal field, particularly for large fl. 

We start by integrating (j3.2.ip over t £ [S', S' -|- T], obtaining 


u{S + T)-u{S) 
T 


QBAs,tu — As,tP^ uU, 


(3.2.3) 


or 



(3.2.4) 
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We want to show that if the left side of (|3.2.4p is small (in some sense) then 
As,tu is close to being a zonal field. 

To do this, we will produce an operator IT with the property that, for 
each s G M, p G (1, oo), IT is a projection of 


= {uG A^) : = 0} onto V^'P n Ker B. (3.2.5) 

As in ()1.1.7p . 

V^’P n Ker .B = {*df : f G H^+^’P{M), Zf = 0}, (3.2.6) 

where, as in (ll.l.dp . 

Z = JVy. (3.2.7) 

As noted in (|3.1.5I) . the geometrical hypothesis on M implies Z = ‘hXs for 
some nowhere vanishing $ G C°°{M). Consequently, 


V^'Pf^KevB = {*df : / G X^f = 0}. (3.2.8) 


We will define IT by 

U{*df) = *dUf, 

where 

n : i7"+i’P(M) ^ 77*+i’P(M) 
is the projection onto KerX3 given by 


with 



The following result is the key to exploiting ()3.2.4p . 


(3.2.9) 

(3.2.10) 


(3.2.11) 


(3.2.12) 


Proposition 3.2.1 If M is a body of rotation about the x^-axis, with posi¬ 
tive Gauss curvature, then, for q G (l,oo), s G M, 

11(7 — n)n||/^s-2,q < Cq^s\\Bu\\Ho^i- (3.2.13) 

Given B{*df) = *dAQ^Zf, from (|1.1.6p . and given (13.2.91) . it suffices to 
prove the following. 
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Proposition 3.2.2 In the setting of Proposition 3.2.1, for p S (l,oo), a G 

||(I-n)/||^^.,P <C||Z/||^^.,P. (3.2.14) 

Proof. Given (13.2.71) . it suffices to show that 

11(1 - n)/||i^..p < CIIXs/IIh-.p. (3.2.15) 

The formula (j3.2.11l) implies that 11 is bounded on L^[M) for all p G (1, oo) 
and commutes with Aq. Also A3 commutes with Aq, so it suffices to estab¬ 
lish 

\\{I-n)f\\Lv<C\\X^f\\Lr. (3.2.16) 

To get this, it suffices to construct a bounded map T on such that 


TA 3 = I - n. 


(3.2.17) 


We construct T in the form 

Tg = ^J ip{s)e^^^gds, 

where e^^^g{x) = g{Rsx). We have 

TXsf = ^f Hs)e^''^Xsfds 


1 


fj'{s)e^^^fds, 


provided ifi—Tr) = To get (I3.2.17p . we want 

— iffs) = 27r(i(s) — 1. 

This is achieved by 

'i/’(s) = S -|- TT, —TT < S < 0, 

S — TT, 0 < S < TT. 


(3.2.18) 




(3.2.20) 


(3.2.21) 


Since if G L^(—vr, vr), T is bounded on each LP{M). This establishes ()3.2.17p . 
hence (13.2.161) . hence (|3.2.14p . hence (|3.2.13l) . The proof of Proposition [3W1] 
is complete. □ 


21 


















Applying Proposition 13.2. ll to ()3.2.4D . we have 


II (/ - 


< 




\n\ 


-|r-i||n(5 + T) - n(5)||H»,, + ^5,T||V„n||H.,,}, 

when u solves (|3.2.1I) . Now (cf. [12], Chapter 17, (2.23)), 
div n = 0 =► = div(n ® u), 


so 


l|V„(t)'u(t)||//-,-j < C'||n(t) (g) n(t)||j:^s+i,,j. 
Meanwhile, as seen in 112.31 with w{t) = rot u{t), 

||u(t)||(,i.oo < C\\w{t)\\L°° < C'(||r(;(0)||Loo + 2|n|), 
hit)\\L2 = ||n(0)||£,2. 

We produce further estimates by interpolation. For starters, 

< c||u(o)|il,"(||w(o)||i« +2|n|)>/2. 

In formulas below, in order to simplify the notation, we set 
O0 = ||n(O)||}^^''(||n;(O)|Ucc+2|O|)^ 
Then (I3.2.26P becomes 

II^(^)IIh1/2,4 < CQi/2- 


Since (|3.2.24p is quadratic in u{t), we want to interpolate ()3.2. 
estimate in (13.2.251) . We get, for 0 < 0 < 1, 


\u{t)\\jjg/ 2 ,pie) < C'||n(t)||^2l'w(i)||^i/2,4 


with 

Now, for 0 < 0 < 1 
with 

r(0) = 


< CQqI2, 


1-9 ^ 9 _2-9 
~~2 ^ 4 ~ 4 ■ 


p{9) 

77®/2’P(®)(M) c 
2p{9) 4 4 


2-9p{9)l2 ^-9 
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2-29 1 - 


(3.2.22) 

(3.2.23) 

(3.2.24) 

(3.2.25) 

(3.2.26) 

(3.2.27) 

(3.2.28) 
with the 

(3.2.29) 

(3.2.30) 

(3.2.31) 

(3.2.32) 























so 


\\u{t)\\j^ 2 /(l-e) < CgLlgi 2 

, 0 < 6» < 1. 

(3.2.33) 

Hence 

||u(t) (g)n(t)||^i/(i-s) 


(3.2.34) 

so, by (3.2.24), 



(3.2.35) 


This indicates taking s = — 1 in (j.S.2.22p . and leads to the following result. 


Proposition 3.2.3 Let u be a global solution to (ll.0.3p . with initial data in 
s > 2. Then there exists Cg < oo, independent of S,T,Ll G (0,oo), 
such that 

11(1 - U)As^Tu\\H-3Me) 

< + r)- «(S)IIh-.-,(.. (3.2.36) 

+ CM0)|li/{|K0)lli»+2|fi|)‘'}, 
for 0 < 6 < 1, with q{9) as in ()3.2.35p . 

Remark. We have 

MS + T) - n(S)||^-i„(,) < Cg\\u{S + T) - uiS)M 

<2Cg\\u{0)\\L2. 

Let us look at some special cases to which Proposition 13.2.31 applies. 
First, if it is a zonal field, then, as seen in 113.11 it is a stationary solution 
to ()1.0.3I) . and consequently the left side of (I3.2.36P vanishes. By contrast, 
recall the non-zonal stationary solutions on S'^ given by (13.1.lip , i.e., 

u = --^X3 + JVgk, (3.2.37) 

Ak — 2 

with Afc = It^+k > 2 an eigenvalue of —A and gk a non-zonal Afc-eigenfunction 
on as in (I3.1.10p . In such a case, 

As,tu^^^ ~ = *dgk. (3.2.38) 

To make contact with the estimate (j3.2.36p . let us suppose that \\u\\l^ = 
||it(0)||^2 Ri 1. Then 

—and ||A^/"(7fc||L2 <(7. (3.2.39) 

A2 — ^ 


23 























It follows that 


\\{I — U)As,TmH-‘‘ ^ C\\dgk\\H-= 


(3.2.40) 


For s € (0,1), this estimate is stronger than ()3.2.36F but of a similar flavor. 
Of course, since (|3.2.37p covers only a special class of stationary solutions to 
(ll.O.ip . it is not surprising that estimates here are better than the general 
estimates guaranteed by (I3.2.36p . 


3.3 Another conservation lavv^ 

As usual, M C is a surface, diffeomorphic to with positive Gauss 
curvature, and invariant under the group of rotations about the xs-axis 
generated by A 3 . As a consequence. 


X is a smooth function of X 3 and 


-— > a > 0 for X G [— 1 , ll. 
dX3 


(3.3.1) 


Next, since A 3 generates a flow by isometries on M, we have div A 3 = 0 on 
M, so there exists ^ G C°°{M) such that 


JV? = -A 3 . 


(3.3.2) 


Clearly As^ = 0. As a further consequence of our geometric hypothesis, 

^ is a smooth function of X 3 and > a > 0 for X 3 G [—1,1]. (3.3.3) 

dX3 

(If M = then X = ^ = X3.) 

We now aim to establish the following conservation law. 


Proposition 3.3.1 Under the hypotheses on M made above, if u{t) solves 
(fFOAD and rot u = w, then 

J f,{x)w{t,x) dS{x) is independent of t. (3.3.4) 

M 
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Proof. From the vorticity equation (I2.2.5p . we have 


fL 

dt 


j dS = j 

M M 

= J ui^X - w) dS (3 3 5^ 

I ^{Vuw)dS + n I fV^xdS. 


M 


Note that (I3.3.ip - ()3.3.3p imply ^ is a smooth function of y; write ^ = ^(x). 
Then fVuX = '^uG{x) where G'{x) = f,{x)- Hence 

J fVuX dS = j VuGix) dS = 0, (3.3.6) 

M M 

since divw = 0 implies V„ is skew adjoint, and V„1 = 0. Next, 

- j Ci^uw) dS = j{VuOwdS 

M M 

= J{JVf,Vf){Af)dS 

M 

= J\x3f){Af)dS 

M 

= {Xsf,Af). 

Now, since X^ commutes with A and is skew-adjoint. 


It follows that 


proving Proposition 13.3.11 


dt 


(3.3.7) 


(A 3 (-A) 1 / 2 /, (-A)V 2 /)= 0 . 

( 3 . 3 . 8 ) 

o' 

II 

( 3 . 3 . 9 ) 

M 



□ 
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3.4 Computation of x and ^ 

Let the surface of revolution M C be given by 



xl + xl = r{x3f, 

(3.4.1) 

i.e., u{x) 

= 0 with u{x) = x\ + x"^ — r(x 3 )^. We have 



Vtt(x) = 2 (xi,X 2, -r(x 3 )r'(x 3 )), 

(3.4.2) 

so the unit outward normal to M is 


N{x] 

Vtt(x) 1 / N // 

' = IV7 ^ M = . ^ n -r(x 3 )r(x 3 )). 

Vu(x)| r(x 3 )y'l + r'(x3)2 

(3.4.3) 

Hence 


(3.4.4) 


We next look for ^ G C°°{M), satisfying = —X^. Clearly ^ is to 
be a function of X 3 , and then the desired condition is |V^| = IX 3 I = r(x 3 ). 
Now 

Ve = ^'(X3)VX3, (3.4.5) 

and Vx 3 is the orthogonal projection onto T^M of 63 , so 

|Vx3|2 = 1 - (63 • N{x))^ = 1 - X{X3)^- (3.4.6) 

Hence ^ is defined by the condition 


^'(^ 3 ) 


rjxs) 

V^-Xixs)'^ 


Bringing in (|3.4.4I) . we obtain 


hence 


so 


Xixsf = 


r'jx^f 
1 + r'{x3Y' 


^-Xix^f = 


1 + r'(x3)2 ’ 


^'(xs) = r(x3)A/l + r'(x3)2. 


(3.4.7) 


(3.4.8) 


(3.4.9) 

(3.4.10) 
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Note that this yields an interesting geometrical interpretation of Namely, 
up to an additive constant, 27r^(x3) is the area of 

{{x,y,z) £ M : z < X 3 }. 


Special case: ellipsoids of revolution 


We specialize our calculations to the case where M is given by 



x^ + X 2 +^—^ = 1 , a > 0 . 

(3.4.11) 

Thus r(x 3 ) = a/1 — 

xi/a? in (j3.4.ip. It follows that 



/(.3) = -^|(i-4)"''. 

(3.4.12) 

hence 

l + r'(x3)2- 
i + r(X3j 

(3.4.13) 

with 

1 1 

a"' a* 

(3.4.14) 

Thus, bv (13.4.41). 

Xixs) = 2 n - 

^1 - /3xi 

(3.4.15) 

and. bv (13.4.lUp. 

^{x3) = 1-15x1. 

(3.4.16) 

For these ellipsoids, X 3 £ [—a,a], and we have < 1, so 

(13.4.15h-(13.4.16h clearlv exhibit y and f as elements of C°°([- 

the formulas 
-a, a]). Note 


that 


0<a<l^/3<0, a=l^/3 = 0, a > 1 ^/S £ {0,1). (3.4.17) 

In case M is the unit sphere, so a = 1, we get 

X{x3) = ^{xs) = xs, (3.4.18) 

as expected. Ellipsoidal planets that bulge at the equator have a < 1. 


27 

















3.5 Smoothness issues 


As we have seen, (I3.4.15p - (|3.4.16l) exhibit x and ^ as elements of a, a]) 

when M is an ellipsoid of the form (13. 4. lip . To extend this, assume 

X3 : M — \—a, al, a = max X3, —a = min 2:3. ( 3 . 5 . 1 ) 

M M 

We claim that if M C is a surface of revolution about the xs-axis, 
diffeomorphic to 5^, and with positive Gauss curvature, the following holds: 

If ‘h E C°°(M) is invariant under rotation about the xs-axis, then 

(3.5.2) 

<l)(x) = (/?(x 3 ) with G C'°°([—a, a]). 

Recall that X )'6 £ C°°{M). We first show how, under the condition (13.5.21) . 
the conclusion x-,^ ^ C'°°([—a, a]) is manifested in the formulas (13.4.41) and 
(j3.4.10p . First note that, by (|3.4.ip . belongs to C^{M), so, under 

the condition (j3.5.2p . 

rixsf = p{x3), p e C°°{[-a,a]). (3.5.3) 

Let us bring in the hypothesis that the curvature of M is nonzero at the 
poles ( 0 , 0 , ±a), so 

p'{±a) / 0. (3.5.4) 

Note how these results can be directly verified in case (13.4.111) . Generally, 


we have 

r’{x3) = ^p'{x3)pix3) 

(3.5.5) 

Hence 


(3.5.6) 

Thus, bv (I3.4.4p. 

V'4/)(x3) + p'(x3)2 ’ 

(3.5.7) 

and. bv (13.4.101). 

= ^V4p(x3) +p'(x3)2. 

(3.5.8) 


By virtue of (I3.5.4p . these formulas clearly give Xi^' G C°°{[—a,a]). 

Geometric hypotheses guaranteeing that the condition (I3.5.2P holds are 
pretty straightforward away from the extreme values X 3 = ±a. Let us verify 
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(j3.5.3p under the following explicit hypothesis on M near the poles (0,0, =ba). 
Namely, we asume M is given near the poles as 


X 3 = ±a + ‘p±{x‘l + X 2 ), 

(3.5.9) 

for Xi+X 2 < 6, with 

99±gC°°([0,<5)), 9 ^±( 0 ) = 0 , ^99'±(0)>0. 

(3.5.10) 

Then, bv (13.4.11). 

X 3 = ±a + ip±{p{x 3 )), 

(3.5.11) 

so 

1 = T±{p{x3))p'{x3). 

(3.5.12) 

Hence 

dp 

(3.5.13) 


which yields p satisfying (I3.5.3l) - (l3.5.4p near X 3 = ±a. 

After these observations, we are now ready to prove a clean smoothness 
result. 


Proposition 3.5.1 Let M C be a smooth, compact surface, invariant 
under rotation about the x^-axis. Assume (I3.5.ip holds. Assume that the 
Gauss curvature K{x) > 0 for all x G M. Then (|3.5.2p holds. That is, 
i/G C°°{M) is invariant under rotation about the x^-axis, then <l>(x) = 
^p{xs) with ip G a,a]). 

Proof. The conclusion about <l>(x) = (fixs) is straightforward except for 
smoothness at X 3 = ±a, so we concentrate on that. Near the poles (0,0, ±a), 
(xi,X 2 ) serves as a smooth coordinate system on M, so 

^{x) = ip±{xi,X 2 ), (3.5.14) 

with 'iIj± smooth on a disk Ds{0) C and invariant under rotations. It is 
a very special case of results of [Ma] that 'ip± are smooth in xf + , so 

$(x) = 7 ±(x?+ xi), 7 ± e (3.5.15) 

This observation applies in particular to X 3 G so we have (I3.5.9p . 

with if^. as in (13.5.101) . The last item of (I3.5.10p . =FV^±(0) > 0, follows from 
X > 0 at the poles of M. Thus the analysis (I3.5.1ip - (l3.5.13p applies, and 
we get 

p{x 3 ) = xl + xllj^ ^ p e C°°{[-a, a]), p'{±a) / 0, (3.5.16) 
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so, near the poles ( 0 , 0 , =ba) 


(3.5.17) 

□ 


$(a;) = -f±{p{x3)), 


smooth in X3 G [—a,a]. 

We take a further look at an ingredient in the proof of Proposition 13.5.11 
namely the following. Let ^^(O) = {{xi^X 2 ) G < (5^}. 

Lemma 3.5.2 // 1 /; G C'°°(Zl 5 ( 0 )) is invariant under rotations, then there 
exists 7 G C'°°([0, (5^)) such that 

V’(xi,X 2 ) = 7 (xi+xi). (3.5.18) 

We present a direct proof of this, not appealing to the general (and 
rather deep) work of [Ma]. It is clear that, if V’ is rotationally invariant, 
then (|3.5.18p holds with 

7 ( 5 ) = 0), sG[0, 5^). (3.5.19) 

The crux of the matter is to show that such 7 is C°° on [0,5^), and of course 
such smoothness is clear except at s = 0. To restate (13.5.191) . we have 

7 ( 5 ) = ■0(si/2), with ^(t) ='i/’(L 0)- (3.5.20) 

We have 

^i-t) = m, (3.5.21) 

and we want to deduce from this that 7 is smooth at s = 0 . 

Now (I3.5.2ip implies that the formal power series of ijj has the form 

00 

(3.5.22) 

k=0 

with only even powers of t appearing. Consider the formal power series 

00 

(3.5.23) 

k=0 

A theorem of Borel guarantees that there exists 7 G C'°°((—<5^, <5^)) whose 
formal power series is given by (|3.5.23l) . Thus 7 (t^) = 'i/’(t, 0) and 7 (t^) both 
have the same formal power series, namely ()3.5.22p . Thus 

l{t^)-l{t^)=u{t), (3.5.24) 
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with 


(3.5.25) 


u e C°^{{-6,6)), u^^\0) = 0, Vj. 

It then follows from the chain rule that 

n(s) = ^ ^ g c 7 ooqq^ (3.5.26) 

Since 

7 (s) = 7 (s)+ u(s), (3.5.27) 

this proves the desired smoothness of 7 at s = 0 . 

4 Stability of stationary solutions 

In this section we examine stability of stationary zonal solutions of (ll.U.ip . 
again assuming M is radially symmetric and has positive Gauss curvature. 
First, in 114. 11 we look at an Arnold-type approach to stability, bringing in 
functionals 

T-l-iu) = j+ ^p{w - fly) dS, (4.0.1) 

M 

for various functions and real constants 7 . Given a stationary solution 
JV/ and w = A/, we see that if 

w{^) — nx(0 is strictly monotone in (4.0.2) 

then one can find and 7 such that u = JV/ is a critical point of (j4.0.ip . 
with positive dehnite second derivative. Stability in H^{M) is a conse¬ 
quence. Note that, for hxed / (hence fixed w), (|4.0.2p holds for all suffi¬ 
ciently large II. 

In ^4.21 we linearize (jl.O.ip about a stationary zonal solution JV/. More 
precisely, we linearize the associated vorticity equation, obtaining a linear 
equation of the form 

^ = rc, rc = -Vjv/C + Vjv(.-qx)^”'C- (4.0.3) 

The symmetry hypothesis on M allows us to write 

r = 0rfc, Tk-.Vk^Vk, Vk = {C€L\M):XsC = ikC}, (4.0.4) 

k 

and deduce that F has spectrum off the imaginary axis if and only if some 
Ffc (A: 7 ^ 0) has an eigenvalue off the imaginary axis. In this setting, we 
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derive a version of the Rayleigh criterion, namely, if T has an eigenvalue 
with nonzero real part, then 

w'{^) — ^x'iO must change sign. (4.0.5) 

Note how this interfaces with the criterion (14.0.21) for Arnold-type stability. 
We see that the Arnold-type criterion for proving stability and the Rayleigh- 
type condition for the lack of proof of linear instability are almost equivalent. 

This is not at all to say that the criterion (14.0.21) nails stability. Just 
when stability holds and when it fails remains a subtle question. The rest 
of this paper is aimed at formulating some attacks on this queston. In (14.31 
we set things up for some specific calculations, which will be continued in 
^ At this point, we will want to make use of classical results on spherical 
harmonics, so in ^4.31 and ^we will specialize to the case M = S"^. 

In ^4.31 we look at (j4.0.4p with 

Tk = ikMk, = MC = A{x^)C + B{x^)/\-\. (4.0.6) 

In the setting of (14.21 ^(xa) = f'{x^) and B{x 2 ,) = — w'{x^). We present 

some results on Spec M^, particularly when 

A{x 3) = af'ixs), B{x 3) = Q. + X^afixs). (4.0.7) 

These results will have further use in ^ 

4.1 Arnold-type stability results 

We use the following variant of the Arnold stability method (cf. [T], pp. 89- 
94, [9], pp. 106-111) for producing stable, stationary solutions to the 2D 
Euler equations, in case M is rotationally symmetric, and has positive Gauss 
curvature. Namely, we look for stable critical points of a functional 

n{u) = + ip{w-nx) + X^w'^dS, (4.1.1) 

M 

with w = rotu and if and 7 tuned to the specific steady solution u. The 
functions y and ^ are as in (11.0.21) and (13.3.21) . See also (I3.4.4P and (13.4.101) . 
Such a functional is independent of t when applied to a solution u{t) to 
(jl.O.ll) . Taking 

u = JV/, so w = A/, (4.1.2) 
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we rewrite (|4.1.ip as 

Hif) = |{i|V/|2 + (^(A/-Ox)+7CA/}d5. (4.1.3) 

M 

Then 


so 


dsH{f + sg) = j |(V/, V5) + s|V5p + y?'(A/ + sAg - ^x )^9 

M 

+ 7 ?Afi(| dS, 

(4.1.4) 


dsH{f + sg)\^^^ 

= j[{Vf,Vg)+ip'{M-^x)^9 + li^9}dS 

M 

= j[-M + -^x) + i^i}9dS. 

M 


(4.1.5) 


This is 0 for all g if and only if / — (p'{Af — fix) ~ 7 ^ is constant, and since 
the stream function / is determined only up to an additive constant, we can 
write 


/ = ¥5'(A/ - fix) + 7 ^, 


(4.1.6) 


as the condition for / to be a critical point of H in (|4.1.3p . Note that (14.1.61) 
implies that, if 

V/IIVC, (4.1.7) 

then 

V(A/- fix) II V/, (4.1.8) 

hence 

(JV/,V(u;-flx)) = 0, (4.1.9) 


so by Proposition 13.1.11 such / produces a stationary solution to (jl.O.ip . 
provided / is a zonal function. If / is not a zonal function, one would 
need to take 7 = 0 in (14.1.11) in order for (j4.1.9p to hold. (Consequently, 
the Arnold method apparently produces much weaker stability results for 
non-zonal stationary solutions than for zonal stationary solutions.) 

To proceed, we apply dg to (|4.1.4I) and evaluate at s = 0, to get 


d^gHif + sg)l^^ = J{\Vg\^ + (^"(A/ - Ox)(A5)2} dS. (4.1.10) 

M 
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Now, if we are given a zonal function /, we want to find ip such that 
(j4.1.6p holds, and then check ()4.1.10l) to see if this is a coercive quadratic 
form in g. Let us write also 


^f = w{0, x = x{0- 

Then (14.1.61) takes the form 

/(O = - ^x(O) + 7^, 


or 


(4.1.11) 

(4.1.12) 

(4.1.13) 


- nxiO) = /(O - 76 

Given arbitrary 7 G M, this identity uniquely specifies ip', provided 

w{^) — nx(0 is strictly monotone in (4.1.14) 

that is, 

w'{^) — is bounded away from 0. (4.1.15) 

With ip' determined, in turn ip is determined, up to an additive constant, 
which would not affect the critical points of (j4.1.3h . Then, applying d/d^ to 
(14.1.131) yields 


<p"{w - nx) = 


Substitution into (I4.1.10p gives 


7 - no 

^x '(0 - ^'( 0 ' 


«fl(/ + .9)U= + (4.1.17) 


(4.1.16) 


M 


f^x'(0 - uj'iO 


By calculations of S ^3. 41(331 as long as the Gauss curvature of M is every¬ 
where positive, both x and | are smooth, strictly monotonic functions of 2 : 3 , 
with positive xs-derivatives, so 

x'(0 > a > 0 on M. (4.1.18) 


As long as the hypothesis (I4.1.14p - (l4.1.15p holds, then either 

^x'{^)-w'{0>b>0, or 

f^x'(0-«^'(0 <-fe<o, 


on M. In the hrst case, we can make 


m) 


7 - r(o 

^ x '{0 - w '(0 


>00 


(4.1.19) 


(4.1.20) 
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on M by taking 7 > 0 large enough, and in the second case we can arrange 
(j4.1.20l) by taking 7 sufficiently negative. Both cases yield 


d^,H{f + sg)l^^>\\Vgf + C\\Ag\\l„ (4.1.21) 

with C > 0, for all 5 G This implies stability of / in as 

a critical point of (14.1.31) (recall that / is defined only up to an additive 
constant), hence stability of u in H^{M) as a critical point of (14.1.111 . We 
summarize. 

Theorem 4.1.1 Given a smooth /(^), u = f is a stable stationary so¬ 
lution to (1-1), in as long as 11 is such that (14.1.14|) - (|4.1.151) hold, 

where w = A/. 

Note that w = A/ implies 

u;(0=/'(0Ae + /"( 6 |Ve|^ (4.1.22) 


if/ = /(0- 

4.2 Linearization about a stationary solution 

Let M C be a compact surface, rotationally symmetric about the xs-axis, 
with positive Gauss curvature, and let u = JV/ be a stationary solution to 
(jl.O.ip . We derive an equation for the linearization at u. More precisely, we 
work with the vorticity equation (I2.2.15p . i.e., 

^ + (JV/,V(u;-L!x)) = 0. (4.2.1) 

Let us set 

fe{t) = f-^sv{t)-\ -, We{t) = w£C{t)-\ -, C = Ar/. (4.2.2) 

Inserting these into the analogue of (I4.2.ip . using (j4.2.1|) and discarding 
higher powers of e produces the linearized equation 

dtC + {JVf, VC) + (JVr?, V{w - Qx)) = 0. (4.2.3) 

Now 

{JVg, V{w - nx)) = -{Vri, JV{w - lly)) ^ 

= .j\j{w-nx)h- 
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Also, since C integrates to 0 on M, we can write 


7? = A-^C, (4.2.5) 

where, here and below, we define to annihilate constants and to have 
range orthogonal to constants. Then (14.2.31) becomes the linear equation 

I = rc. (4.2.6) 

where 

rC =-Vjv/C + Vjv(u,-Qx)^~^C. (4.2.7) 

The question of linear stability is the question of whether T generates a 

uniformly bounded group of operators on 

Ll{M) = {c e : J CdS = 0 y (4.2.8) 

M 

Under our hypotheses, we have x = x(0) with ^ as in (13.3.21) . i.e., 
JV^ = —X 3 . Let us also assume / is a zonal function, i.e., X^f = 0, so 
/ = This also implies X^w = 0, hence w = w{^). Then 

JV/ = -/'(02f3, 

JV(7/;-L!x) = [L!x'(e)-«;'(e)]2^3, 

and (14.2.71) becomes 

rc = /'(O^sC + (^2x'(6 - w;'(0)^3A-ic. 


In such a case, T commutes with A 3 , hence we can decompose 


Ll{M) = ^Vk, 

k 

(4.2.11) 

where, for A: G Z, 


Ufc = {C G lUm) : AgC = ikC}, 

(4.2.12) 

and we have 


r = 0rfc, Tk-.Vk^Vk, 
k 

(4.2.13) 

where 


TfcC = ik[fm + (f2x'(0 - 

(4.2.14) 


(4.2.9) 

(4.2.10) 
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Note that 


A ^ : Vfc —> Vfc is compact 


(4.2.15) 


for each k, so each is a compact perturbation of a bounded, skew-adjoint 
operator on 14. In light of this, basic analytic Fredholm theory yields the 
following. 

Proposition 4.2.1 For each k, 


(4.2.16) 


SpecFfc C ikT; U Sk, 


where 


S = {f'{X) : oo < < cki}) «o = ™in C; (4.2.17) 


and Sk is a countable set of points in C whose accumulation points all must 
lie in ikTi. Each p, G Sk is an eigenvalue ofTk, and the associated generalized 
eigenspace is finite dimensional. 

In fact, for each p £ C \ ikTi, Tk — pi is a bounded operator on 14 that 
is Fredholm of index 0, and it is clearly invertible for \p\ > ||rfc||. 

Corollary 4.2.2 Assume F has the form (I4.2.10jl . If SpecF is not con¬ 
tained in the imaginary axis, then some F^ has an eigenvalue with nonzero 
real part. 

Now having Spec F C would not guarantee that F generates a bounded 
group of operators on L^(M), but not having this inclusion dehnitely guar¬ 
antees that the associated group of operators is not uniformly bounded. 
Thus Corollary 14.2.21 points to an approach to finding cases that are linearly 


unstable. 


Actually establishing such cases of linear instability is not so straight¬ 
forward. We proceed to derive some necessary conditions for such linear 
instability to hold, i.e., for some F^ to have an eigenvalue with nonzero real 
part. 

Of course, Fq = 0. Suppose /c / 0 and F^ has an eigenvalue p = ikft, /? ^ 
M. Then there exists a nonzero C € 14 such that 


(/'(O - /3)C = -(f^x'(0 - 


(4.2.18) 


hence 



(4.2.19) 
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where r] = Note that if /3 ^ M, the denominator on the right side of 

(14.2.191) is nowhere vanishing. In (|4.2.18l) - (l4.2.19p . Q and r/ would not be 
real valued. Taking the inner product of both sides of (14.2.19p with r] yields 


{Ar],r]) 



w’ic) - nx'ic) 


\ri?dS 



w'{0 - 

l/'(0-/3P 




(4.2.20) 


Now (Ary, ry) is real and negative, but Im/3 ^ 0. Hence taking the imaginary 
part of (I4.2.20P yields 


fw'iO-Qx'iO 

J l/'(0-/3P 


|ryp dS = 0. 


(4.2.21) 


Using this in (|4.2.20ll gives 

(Ary,ry) = J [/(O - K] |ryp dS < 0, (4.2.22) 

S2 


We have (14.2.211) and (I4.2.22P as necessary conditions for T^ to have an 
eigenvalue with nonzero real part, with associated eigenfunction ( = Ary, ry G 
14 . These results in turn imply the following. 


Proposition 4.2.3 IfT has an eigenvalue with nonzero real part, then 

w'{s) — Q,x'{s) must change sign in s G (aojCti)) (4.2.23) 

with Oj as in (14.2.171) . and 

VK € M, 3sG(aO)«i) such that {w'{s) — nx'is)){f'{s) — K) < 0. 

(4.2.24) 

In the setting of planar flows (and with H = 0), (14.2.231) is known as the 
“Rayleigh criterion” for linear instability, and ^4.2.241) is called the “Fjortoft 
criterion.” See [9], pp. 122”123. 

Proposition 14.2.31 is close to Theorem 14.1.11 in the following sense. By 
Theorem 14.1.11 if 

rr;^(s) — Hx^('S) / 0 for all s £[ 00501 ]) (4.2.25) 

then the associated stationary solution u = JV/ to (jl.O.ip is stable, in the 
sense of Condition (14.2.231) is a little stronger than the assertion that 

(j4.2.25p fails. Thus, in some sense, the hrst part of Proposition 14.2.31 is 
almost a corollary of Theorem 14.1.11 
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4.3 Further results on linearization 

Here we produce some results complementary to those of 1I4.2I We consider 
operators F of a more general nature than those in 1)4.21 as indicated in 
(I4.3.4P below. However, we specialize from more general surfaces of rotation 
to the standard sphere S'^, in order to make some explicit computations 
using spherical harmonics. 

To proceed, we investigate matters related to whether the operator T 
generates a uniformly bounded group on = {/ G Lp‘{S‘^) : Jg 2 f dS = 

0 }, when T has the following structure: 

r = 0rfc, Fk-.Vk^Vk, Vk = {f£Ll{S^):Xsf = ikf}, (4.3.1) 

k 

where is the vector field generating 27r-periodic rotation about the X 3 - 
axis. We assume 

Tk = ikMk, Mk = M\y^, (4.3.2) 

where 

MC = A(x3)C + B{x^)X-^C- (4.3.3) 

We assume A and B are smooth and real valued. In studies of linear stability 
of stationary, zonal Euler flows on the rotating sphere, such an operator 
arises with 

MX 3 ) = f'ixs), B{x 3 ) = n-w'{x 3 ), (4.3.4) 

with w = Af = rot u, u a steady zonal solution to the Euler equation. 

The question we examine is whether Spec T is contained in the imaginary 
axis. In view of (14.3.21) . this is equivalent to the question of whether SpecM^ 
is contained in the real axis, for each k ^ 0. Basic Eredholm theory gives 
the following. (Compare Proposition 14.2.11 1 

Proposition 4.3.1 For each k ^ 0, 

SpecMfc C S U 

where S = {^(xs) : —1 < X 3 < 1} and is a countable set of points in C 
whose accumulation points all must lie in S. Each X € is an eigenvalue 
of Mk, and the associated generalized eigenspace is finite dimensional. 

In fact, for each A G C \ S, Mk — A/ is a bounded operator on Vk that is 
Eredholm, of index 0, and it is clearly invertible for |A| > ||Mfc||. The next 
result is a cousin to Corollary 14.2.21 
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Corollary 4.3.2 Assume M has the form (j4.3.3l) . If SpecT is not con¬ 
tained in the imaginary axis, then some {k ^ 0) has an eigenvalue that 
is not real. 


Remark. If A is a non-real eigenvalue of then A is an eigenvalue of 
both Mk and 


Actually establishing such cases of linear instability is not so straight¬ 
forward. We proceed to derive some necessary conditions for such linear 
instability to hold, i.e., for some to have a non-real eigenvalue. 

If A ^ M is an eigenvalue of , then there is a nonzero Q ^Vk such that 


(A(x3) - A)C = -B{x3)A ^C- 


We can take g 


A £ Vk and write this as 


Ag 


Bjxs) 

A - A{x3)^' 


(4.3.5) 


(4.3.6) 


Note that if A ^ M, then, since A is real valued, the denominator on the right 
side of (14.3.61) is nowhere vanishing. (Note also that g is not real valued.) 
We take the inner product of both sides of (I4.3.6P with g, to get 




Now (A? 7 , g) is real and negative. Hence the imaginary part of the last 
integral is zero. If A ^ M, this forces 




Given this, we can then deduce from (4.3.6A) that 

i^V,v)= J j^^^j^{K-A{x3))\gfdS<0, VReM. (4.3.9) 
52 

We have (I4.3.8P and (|4.3.9I) as necessary conditions for to have an 
eigenvalue A ^ M, with associated eigenfunction C, = Ag, g £ Vk- These 
results imply the following. (Compare Proposition 14.2.31 1 
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Proposition 4.3.3 IfT has a non-imaginary eigenvalue, then 

B{s) must change sign in s £ (—1,1), (4.3.10) 

and 

V K S M, 3 s S (—1,1) such that B{s){K — ^(s)) < 0. (4.3.11) 

Condition ()4.3.10p is a version of the “Rayleigh criterion” and (14. 3. lip a 
version of the “Fjortoft criterion” for linear instability. Compare the remarks 
after Proposition 14.2.31 

Regarding the relation between (14.3.101) and (|4.3.1ip . we mention that 
there is at least one situation where the Rayleigh criterion (14.3.lOp holds but 
the Fjortoft condition (|4.3.11l) fails, namely when A{x^) = A is constant. 
Then (14.3.lip fails for K = A, but (14.3.lOp holds for many choices of B^x^). 
This result is equivalent to the statement that 

B{x 3 )A~^ has real spectrum on each 14, (4.3.12) 

for k ^ 0. This fact might seem nontrivial, since R(a: 3 )A~^ is not self adjoint 
(if B{x‘i) is not constant), but this operator acts on Sobolev scales, and in 
this framework the operator is similar to, and has the same spectrum as 

- (-A)-^/2^(x3)(-A)-i/2^ (4.3.13) 


which is self adjoint. 

Regarding the reverse implication, we have: 

Proposition 4.3.4 Assume A and B are continuous and real valued on 
[-1,1]. Then (14.3.111) ^ (14.3.101) . 

Proof. There exist Ki and K 2 such that Ki — A{s) > 0 for all s £ [—1,1] 
and K 2 — 4(s) < 0 for all s £ [—1,1]. Applying (|4.3.1ip to K = Ki yields 

51 £ (—1,1) such that B{si) < 0 and applying (14. 3. lip to iF = K 2 yields 

5 2 £ (—1,1) such that B{s 2 ) >0. □ 

Here is another case where (|4.3.10p does not imply (|4.3.1ip . Namely, 
A(s) = —B{s), where (|4.3.11ll fails for K = 0. 

R seems not so easy to give examples where (|4.3.10p holds but (14.3.111) 
fails when A(s) and B(s) have the form (|4.3.4p . with / and w zonal functions 
related by u; = A/. Suppose, for example, that / is a zonal eigenfunction 
of A, 

A/ = -Kf, so w = -\yf {\y > 0). (4.3.14) 
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Then 


B{s){K - A{s)) = {n + Kns)){K - fis)) 

= -X.[-Y^-ns)){K-f{s)). 


(4.3.15) 


To pick K £ U. violating (j4.3.1ip . we need the two factors above to have 
the same zeros, to avoid the product changing sign. This tends to force 
K = —Vl/Xy. But then B{s){K — A{s)) = —X^{K — /'(s))^, which is < 0 on 
most of (—1,1). So this approach fails to produce an example where (j4.3.10p 
holds but (j4.3.1ip fails. 

Having the Fjortoft condition hold along with the Rayleigh condition is 
certainly an acceptable state of affairs, and it is of interest to pursue the 
use of such / as in (j4.3.14|) . We will find it useful to generalize a little, and 
consider the situation 


A/=-At,/, w = -^if (//> 0 ). 


(4.3.16) 


Note that 

Spec(-A) = {Xj = j{j + 1) : j = 0,1, 2, 3,... }, 
Spec(-A)|^^ = {Xj : j > \k\}. 

Let us take 

/(xs) = aPuixs), d > 0 , 
where Piy are Legendre polynomials, given by 

1 




for example. 


Po{s) = l, Pi{s) = s, P 2 {s) =-1), 

Ms) = - 3s), PM) = ^(35s^ - 30^2 + 3). 

Z o 


(4.3.17) 

(4.3.18) 


(4.3.19) 


Taking / = aPo produces a trivial flow. Taking / = dPi gives f'{s) = a, 
hence w'{s) = —QjiOLi constant, so B{s) = Xl — w'{s) does not satisfy (I4.3.10p . 
The hrst choice that might lead to linear instability is 


/(xs) = aP2{x3), 


(4.3.20) 
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giving 


fixs) = ax3 


w'{x3) = -^ax3 


(4.3.21) 


(with a = Sd), hence 

^( 3 ^ 3 ) = 01 x 3 , B{x 3 ) = O + fiax 3 . (4.3.22) 

Then, given 12 >0, the Rayleigh condition ()4.3.10p holds if and only if 

0<n< fia. (4.3.23) 


For n > fia, the Arnold stability criterion applies. Hence, by a limiting 
argument, for H = /ra, F will have no non-imaginary eigenvalues. 

When (j4.3.23p holds, we might find that F does have some non-imaginary 
eigenvalues. That is, might have some non-real eigenvalues, for some 
k 0. Let us take a closer look at this issue when H = 0. In such a case. 


AffcC = ax3 + /raxsA 
= ax3{I + 


(4.3.24) 


for C, ^Vk- Recall we are assuming /i > 0. Now, by (j4.3.17p . 

0 < /r < Afc I -|- /iA“^ is positive definite on I 4 

=> SpecMfc = a Spec(/ -|- /rA“^)^/^X 3 (/ -|- 
=> SpecMfc C M. 

(4.3.25) 

A limiting argument gives the last conclusion for ^ = Xk- We record the 
conclusion. 


Proposition 4.3.5 In case A and B are given by (|4.3.22p and H = 0, 

0 < ^ < Xk SpecMfc C M. (4.3.26) 

Now let us specialize to the case relevant for Euler flow. That is to say, 
we take /r = A 2 = 6 in p4.3.21l) - p4.3.22p : 

A{x 3) = ax3, B{x 3) = H + X2ax3. (4.3.27) 


Corollary 4.3.6 In case A and B are given by (|4.3.27p and H = 0, 

k>2^ Spec Mk C M, (4.3.28) 


and ditto for M-k- 
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Thus, if linear instability arises in this situation, the only possibility is 
that 

Spec Ml is not contained in M, (4.3.29) 

and ditto for M_i. It is therefore of great interest to investigate whether 
(14.3.291) holds. 

Let us extend our considerations to nonzero 12, in the setting of (|4.3.27p . 
Then we have 

(Mfc + Ag ^S2L)C = + -^2 + (12 + \2ax5)A ^ 

= (crx3 + A 2 ^f2)[/ + A 2 A ^](^, 

for C G 14 . Again, we see that has real spectrum if A: > 2, so again our 
search for non-real eigenvalues of is reduced to investigating whether 
(14.3.291) holds. 

Keep in mind that Arnold stability holds for 12 > A 2 a, in this situation. 
Thus we are looking at when (|4.3.29l) holds, given 

0 < 12 < A2a. (4.3.31) 


Note. For the purpose of this analysis, there is no loss of generality in 
taking a = 1. 

We next generalize the setting ()4.3.27p . along the lines of ()4.3.14p . Thus, 
in place of (|4.3.27p . we have 

A{x 3 ) = af'{xs), B{x‘i) = 12 X^af'ixs). (4.3.32) 

Now, in place of (I4.3.30p . we have 

(A4 + X-^ni)C = iaf'ixs) + X-^n)C + {n + X,anx3))A-\ ^ 

= {af{xs) + X-^n)[I + X,A-^]C, 

for C G 14 . This is a composition 

X-^B{x3){I + X,A~^), (4.3.34) 

and this operator has real spectrum as long as either factor, B{x3) or I + 
A,yA“^ is positive, as an operator on I 4 . If 12 is such that B{xk) changes 
sign, the operator still has real spectrum as long as / -|- Ai/A“^ is positive 
on 14 , i.e., as long as Xy < Aj^j. This produces the following variant of 
Proposition 14.3.51 
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Proposition 4.3.7 In case A and B are given by ()4.3.32p . then 


SpecMfc C M (4.3.35) 

provided that either Q + X^af'{x 3 ) does not change sign or Aj^ < A|fc|. 

5 Appendix by Jeremy Marzuola and Michael Tay¬ 
lor: Matrix approach and numerical study of 
linear instability 

In 0we saw that the sufficient condition (j4.0.2l) for stability in H^{M) of 
a stationary zonal solution u = JV/ and the necessary condition (|4.0.5I) 
for the existence of non-imaginary spectrum of the linearized operator T 
in (I4.0.3P are almost perfectly complementary. Nevertheless, as we will see 
here, the spectrum of T might be confined to the imaginary axis even when 
(I4.0.5P fails. Equivalently, the operators : 14 —)• 14 in (|4.0.6I) might all 
have real spectrum, even in cases where (I4.0.5P fails. Here we specialize to 
M = S'^ and make some calculations in cases 

f{x) = cPuixs), 1 ^ = 2,3, 4. (5.0.1) 

The operator takes the form 

MkC = cPUx3)C + {n + X,cPUx3))A-\, CeEfc. (5.0.2) 

The spaces 14 have orthogonal bases 

{P'^'/^pHxs) ■■ i > \k\}, (5.0.3) 

which can be normalized to produce orthonormal bases. Classical identities 
for spherical harmonics lead to representations of Mj. as infinite matrices. 
We carry out these calculations for u = 2 in T5.1I and for z/ = 3 in ^5.21 
For short, we sometimes refer to the matrices associated to in (15.0.21) 
as Pu{Vk) models. 

In 115.1l we use the matrix representation of Mi (for n = 2) to prove that, 
for all II > 0, Ml has only real spectrum. (That Mk has only real spectrum 
for I A:I > 2 in this situation follows from Corollary 14.3.61 ) By contrast, 
the Rayleigh-type condition (j4.0.5p guarantees Mi has only real spectrum 
provided H > A 2 = 6, but it does not apply to II G [0,6). This extra 
constraint on the spectrum of Mi for such small II was first suggested to 
the authors by output from a Matlab program. Having seen the output, we 
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were able to prove that such a constraint holds. We also show that Mi has 
a generalized 0-eigenvector at = 0, giving rise to a weak linear instability. 

In go we work out the infinite matrix representations of Mi and M 2 
(for 1 / = 3), acting on Vi and V 2 , respectively. (In this case, Proposition 
14.3.71 implies that M^ has only real spectrum for \k\ > 3.) The Rayleigh- 
type condition (I4.0.5P guarantees that Mi and M 2 have only real spectrum 
provided n > (4/5)A3 = 48/5. The analysis of Mi and M 2 is more difficult 
than that of Mi in 115.11 At this point, we have numerical results on trunca¬ 
tions of these matrices that indicate linear stability for substantially smaller 
values of n than 48/5. 

These numerical results are discussed in ^5.31 There we take N x N 
matrix truncations of the operators M^, arising in (15.0.21) . for v = 
3,4, k < V. After some discussion about stabilization of the non-real spec¬ 
trum of such matrices for moderately large N, we take N = 400. We use 
Matlab to find the non-real eigenvalues and graph their imaginary parts, as 
functions of 11. These graphs indicate linear stability for somewhat less 
restricted than what the Arnold-type stability analysis of 114. II requires. We 
also see numerical evidence of how stability might not be simply a monotone 
function of 12, for 12 > 0. Taken together with the rigorous results we have 
established through ^5.11 these numerical results suggest much interesting 
work for the future. 

5.1 Matrix analysis for /(x) = cP 2 {x^) 

Here we pursue the question of when (j4.3.29l) holds. We recall the setting. 

A4 = M|^^, Pfc = {/e 4(^2) : A 3 /= zfc/}, (5.1.1) 

and 

MC, = A{xs)C + B{x3)A~^C. 

We take 

A{x3) = X3, B{x3) = Q +X2X3, A2 = 6 , 

and ask the following. 

Question. For what values of 12 does 

A4i have a non-real eigenvalue? (5.1.4) 


(5.1.2) 

(5.1.3) 
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We assume > 0. As we have seen, the “Rayleigh criterion” produces 

0 <n<X 2 ( 5 . 1 . 5 ) 

as a necessary condition for (I5.1.4p to hold. We want to see how close (15.1.51) 
is to being sufficient. In the context of (|5.1.3p . it will turn out to be far from 
sufficient. 

To investigate this, it is convenient to represent Mi as an infinite matrix. 
An orthogonal basis of Vi is given by 

0 = e'^P/(x3), £>1. (5.1.6) 

Here P/ is an associated Legendre function given in (7.12.5) of [8] as 


P/(f) = -(i-f 2 )i/ 2 p;(t). 


We mention that 


— £{£ + l)- 

One has m, p. 201, #10) 

pHt? dt = ^ + 

y_i ^ ^ ^ 21 + 1 {£-l)\ 2£+l ' 


(5.1.7) 

(5.1.8) 

(5.1.9) 


Hence, up to a constant, which we can ignore, an orthonormal basis of Vi is 
given by 


Q = 


2£+l 


2 £{£ + 1 ) 

The operator Mi is given by 


£>l. 


(5.1.10) 


MiC£ = A{x^)Q + Bix3)A-^Q 
= A{x3)Ce - 

Xe 

= X3Ce - -r(^ + ^‘2X3)Ci- 


(5.1.11) 


To proceed, we need to write as a linear combination of {Cj}- To 
do this, we use (7.8.4) and (7.8.2) of [8], 


tp',{t) = pLiit)+mt), 

( 2 ^ + i)p#t) = p;+i(t)-p;_i(t), 


(5.1.12) 
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which combine to give 


tP^{t) = 


1 


2£+l 

Plugging this into (I5.1.1U|) yields 




(5.1.13) 


xsCe = 


(2^ + l)(2£- 1) 


0-1 + 


£{1 + 2 ) 


{2e + l){2£ + 3) 


0 + 1 - 


(5.1.14) 


We use the obvious convention that Co = 0- It is illuminating to write this 
as 


xsCi = cLiCe-i + o^+iO+i) = 




(5.1.15) 


{2£ + l){2£ - 1)' 

Thus the matrix representation of C a^sC on Vi has the 3x3 truncation 
0 02 


= 

Returning to (15.1.111) . we have 


02 0 03 j = j yiTs 

03 0 


f ^ 7175 

0 v^8/35 

v/8/35 0 


(5.1.16) 




MiCi — o^O-i + O£+i0+i —^(ofO-i + Of+iO+i) — -^0 

A£ A£ 


— OfO-l + O£+l0+l — 


A 2 


-1 + O£+l0+l + 


A 2 


Recall that 


(5.1.17) 


(5.1.18) 


Xi — £{£ + 1 ), A 2 — 6 , 

and £ runs over {1, 2, 3,... } in (I5.1.17p . In particular, the 3x3 truncation 

Mf = - R^^'>, (5.1.19) 


of Ml is 



with as 

in (|5.1. 

16P, and 


/ 0 

02 /A 2 

= X 2 

02 /A 1 

0 


V 

03 /A 2 


0 /Ai 


I/A 2 


0 




(5.1.20) 


As it turns out, Matlab programs strongly indicate that Mi has only real 
spectrum, even at 11 = 0. Stimulated by such programs, we have managed 
to verify the results they suggest, and prove the following two propositions. 
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Proposition 5.1.1 Take 12 = 0, and let be given by ()5.1.15p . Xi by 
(IS.l.Sp . Then Mi has only real spectrum on Vi. 

Proof. In this case, 


MiCi = -2 x3Ci, M 1 C 2 = 0, 

MiO = (^1 - for 1 > 3. 


(5.1.21) 


The formula (|5. 1.151) for X 3 Q gives 


Range Mi C V 12 = Span M > 2}. 


(5.1.22) 


Thus any eigenfunction of Mi must lie in Vi 2 . Now 


Ml : V12 V12, Ml = X3{I + X2A-^ 


(5.1.23) 


implies 


has real spectrum, 


(5.1.24) 


by the argument proving Proposition 14.3.51 This proves Proposition 15.1.11 


□ 


Notwithstanding Proposition 15.l.H we do have linear instability at 12 = 
0. In fact, it follows from (j5.1.211) and (j5.1.15p that 


= Cl - ^C2 


(5.1.25) 


so : t € M} is not uniformly bounded on Vi. 

We next extend Proposition 15.1.11 to cover the case 12 > 0. 

Proposition 5.1.2 In the setting of Provosition 13.1.11 (i.e., with A^x^) = 
X3, B{x 3 ) = 12 + A2X3J, take 12 > 0. Then Mi has only real spectrum on Vi. 

Proof. In place of (I5.1.2ip . we have 


Mi(i = —2x3^1 + Ai ^12Ci, 

M 1 C 2 = A2-^12C2, 


(5.1.26) 
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It is convenient to rewrite this as 


(Ml - Xi^n)Ci = -2x3Ci, 

(Afi — II)C2 = (-^2 ~ X^ )f^C2) ^ 27) 

(Ml - X^^nXi = (l - ^)x3Ce + (A7^ - X^X^Q, t > 3. 

Now (15.1.151) plus (j5.1.27|) yields 

Range (Mi — C R 12 = SpanjC^ M > 2}. (5.1.28) 


Hence any eigenfunction of Mi must lie in V 12 . Also, (I5.1.28P implies Mi : 
V12 ^ Ri2, so 

Ml - A 2 : Ri2 ^ Ri2. (5.1.29) 

On the other hand (parallel to (I4.3.33I) - (I4.3.34I) . noting that (14.3.321) holds 
with Xv = A 2 ), 

(Ml - A 2 ^H)C = A 2 ^R(a:3)(/ + A2A-i)C, (5.1.30) 


for C G 1 / 12 , and since / + A 2 A ^ is positive semidefinite on I/ 12 , it follows 


that 

This proves Proposition 15.1.21 


has real spectrum. 


(5.1.31) 

□ 


We conjecture linear stability when H > 0: 

Conjecture. In the current setting, : t G M} is uniformly bounded 

on Vi for each H > 0 . 


5.2 Matrix analysis for /(x) = cP^xX) 

We work with the following modification of the setting of ^5.11 As there, 

Pfc = {/ G Ll{SX : A 3 / = ikf], (5.2.1) 

and we set M^ = MI 14 , with 

MC = A{xXC + B{xXIX-X- (5.2.2) 


As before, 


A(x3) =/'(X3), B{x 3 ) = ^-w'{xX), 


(5.2.3) 
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where w = A/. As before, we take /{x^) to be a zonal eigenfunction of A, 
hence a multiple of for some We saw in ^4.31 that taking z/ = 0 

or z/ = 1 does not work, and in 1)5.11 that taking v = 2 does not work, to 
produce examples of non-real eigenvalues of Mi. Here, we take v = 3. Now 


/(xs) = aPsixs) 
General formulas yield 


A/ = -As/, As = 12. 


15 


Psii) = lienee P^{t) =—{t - - 


2 2 

Hence, in (15.2.21) . we will take 


A(x 3 ) = x|-i, H(xs) = 11 + As(xi - 


(5.2.4) 


(5.2.5) 


(5.2.6) 


In light of Proposition 14.3.71 we are interested in the behavior of Mi and 

Ms. 

We start with an analysis of Mi. We take the orthonormal basis {Q : 
^ > 1} of Vi given by (j5.1.6l) - (|5.1.inp . As noted there 



+ 

II 

-< 

>< 

1 

II 

<1 

(5.2.7) 

and (cf. (5.1.15)) 

xsCe — o-£C£-i + a^+iC^+i, 

(5.2.8) 

with 

/(£+!)(£-!) 

^ Y (2^ + l)(2£-l)- 

(5.2.9) 

In the current setting. 




MiO= (xi-i)o-^(ll + As(xi-i))o, (5.2.10) 


A^ 

so it is useful to note that (|5.2.8I) implies 

+ (a| + + a£^ia£^2Ce+2, 

or equivalently 


with 


X^Ci = hCl-2 + qO + ^£-1-20-1-2, 


— a£a£—i, C£ — aj P fl^+i- 


(5.2.11) 

(5.2.12) 

(5.2.13) 
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In ()5.2.7p - (j5.2.12p . i >1. We use the natural convention that Co 
Putting together (j5.2.10p and (I5.2.1ip yields 


MiQi = biQ-2 + {c£ - + bi+2Ci+2 

— ^ (j}eCl -2 + (cg ~ 5)^^ ^^+2C£+2) 

Our goal is to investigate for what H > 0 does 

Ml have a non-real eigenvalue. 

As we know, the “Rayleigh criterion” produces 

0 < n < -A 3 = A 3 max (x^ — , 

5 |a;3|<l V 5/ 


C-1 = 0. 


(5.2.14) 


(5.2.15) 


(5.2.16) 


as a necessary condition for (15.2.15|) to hold. We make a numerical study of 
(j5.2.14p to indicate how close ()5.2.16p is to being sufficient. 

Numerical experiments, described in ^5.31 indicate that (15.2.15p holds 
for 0 < 0 < 7 with 7 ~ 1. This is a lot smaller than (4/5)A3 = 9.6. 

We move along from Mi to M 2 , i.e., we take k = 2 m (I5.2.1|) . Parallel 
to ()5.1.6p , an orthogonal basis of V 2 is given by 

Ci = e^^^Pi{x3), t > 2. (5.2.17) 


In this case, the associated Legendre function is given in (7.12.5) of [ 8 ] 
as 


Pf{t) = {l-t^)P';{t). (5.2.18) 


We mention that 

AQ = -XeQ, Xe = £{£ + !). (5.2.19) 

We emphasize that here £ >2. In (|5.1.6p - (l5.1.8l) . we had £ > 1. As for the 
norms of these functions, one has m, p. 201 , # 10 ) 



2 {£ + 2 )\ 

2£ + l {£-2)\ 


2{£ + 2){£ + l)£{£-l) 
2£ + l 


(5.2.20) 


Hence, up to an unimportant constant, an orthonormal basis of 1# is given 
by 




2£ + l 


2{£ + 2){£ + l)£{£-l) 


e2*^(l-x^)P;(x3), £> 2 . 


(5.2.21) 
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The operator M 2 acts on this basis as 


MsCr = A{x3)Q + Bix3)A-% 


= A{x3)C£ - 


(5.2.22) 




To proceed, we need to write x^Q as a linear combination of {Cj}- Ihom 
(j5.1.12p and (j5.1.13p . we get, upon applying d/dt, 


'+1 


tp^t) + Hit) = ^^PLiit) + 

P'Jt) = — —p',' 1 it) + -^—Pm it), 
> 21 + 1 2i + l ’ 


hence 


tP'Ht) = 


21 + 1 

Plugging this into p5.2.2ip gives 


P'l-ii^) + 7T7TT-' ^+ 1 '- 


/ -I- 2 

' P'L,it). 


21 + 1 


(5.2.23) 


(5.2.24) 


X 3 Q = 


' ie + 2)ie-2) 

i2i + l)i2£-l) 


Ce-i + 


' (£ + 3)(£-l) 
(2£ + 3)(2£ + l) 


0+1, i > 2. (5.2.25) 


We use the natural convention that Cj = 0 for j < 2. It is convenient to 
write (I5.2.25P as 


^ ^ ^ (^ + 2 )(^-2) 

X 3 O - arO-i + a,+iO+i, a, - y ^ . 

(5.2.26) 

This in turn implies 


xlCe = o^ar-iO-2 + (a| + o|+i)0 + «£+ia£+20+2, 

(5.2.27) 

or equivalently 


XsCi = 00-2 + qO + O+ 2 O+ 2 ) 

(5.2.28) 

with 


u 2,2 

0£ — C£ — cl£^i- 

(5.2.29) 
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Recall that in (j5.2.26l) - (l5.2.29p . (.>2, and we use the convention that Cj = 0 
for j <2. Putting together (I5.2.22P with (I5.2.28P yields 


M2Q = — -^0 + ^£+ 20+2 



- ^ {biCi -2 + (^Ci - 0 + ^£+ 20 + 2 ) 

-yCe, 

(5.2.30) 

for i > 2. 



Our goal 

is to investigate for what 11 > 0 does 



M 2 have a non-real eigenvalue. 

(5.2.31) 

As we know, 

the “Rayleigh criterion” produces 



0 < H < ^A 3 = A 3 max fxo — 

5 h3|<i 1 5/ 

(5.2.32) 


as a necessary condition for (j5.2..Sip to hold. We make a numerical study of 
(j5.2.30l) to indicate how close (|5.2.32p is to being sufficient. 

Numerical experiments, described in ^5.31 indicate that ()5.2.3ip holds 
for 0 < < 7 with 7 ~ 1/2, and that (j5.2.3ip ceases to hold for > 7. 

Again, 1/2 is a lot smaller than (4/5)A 3 . 

In ^5.31 we will also examine such matrices that arise when /(xs) = 
cPiix^)- More precisely, we take 

A(x 3 ) = - ^X 3 , R(x 3 ) = n + A 4 (^x| - ^ 3 : 3 ^, (5.2.33) 

and define 

A/fcO = 0 - ^ + A 4 [xl - ) Ci, (5.2.34) 

on the orthonormal basis {C^} of I 4 , given by (|5.1.6p - (j5.1.10p for A: = 1, by 
(j5.2.2ip for k = 2 and by a comparable strategy for k = 3. We make use of 
matrix formulas parallel to (15.2.141) and (|5.2.30l) in this setting, which are 
given in ^5.31 
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5.3 Numerical study of truncated matrices 

Here we study truncated versions of the matrix operators arising from the 
attack described in ^5.2l on stability of banded structures for the Euler equa¬ 
tions with Coriolis forces on the sphere. Before describing how this is done, 
we mention one result that leads one to believe that eigenvectors of such op¬ 
erators Mfc as described in (I5.U.2I) . or more generally (I4.3.2l) - (l4.3.3p . should 
be expected to be captured fairly accurately by such a truncation. 

Proposition 5.3.1 Given A and B real valued and real analytic on 5^ in 
(j4.3.3p . if ( £ Vk is an eigenvector of M^, with eigenvalue A ^ M, then C, is 
real analytic on 5^. 

Proof. As seen in (|4.3.6I) . rj = satisfies an elliptic partial differential 

equation with analytic coefficients, so real analyticity of rj, hence of ^ = Arj, 
follows. □ 

Given that such an eigenfunction C € 14 of is real analytic, its 
spherical harmonic expansion is rapidly (in fact, exponentially) convergent. 
The truncation of such should then be a high order quasimode of the 
associated truncation of M^. We are not currently in a position to derive 
rigorous conclusions about the spectrum of from numerical results on 
truncations, but we are motivated to take such numerical results as a strong 
indication of how the spectrum of Mk behaves. 

Numerically, we implement in Matlab the eigenvalue solver eigs for an 
N X N block truncation of matrices Mk, such as arise for /c = 1 in (j5.2.inh 
and for k = 2 ui (|5.2.22p . We will refer to these finite matrices as . 

Our first observation is that the non-real spectrum of appears to 
stabilize before N becomes particularly large, and to be set in the upper left 
block of such a matrix. Figured] illustrates this for N xN matrix truncations 
for the model, with 0 = 1/5, and for the ^ 3 ( 14 ) model, with 0 = 2. 

Figure [2] has analogous illustrations for the 14(14) model, with O = 1/10, 
and the 14(14) model, with 0 = 1. In three of these four cases, one sees 
stabilization of the non real spectra in truncations well before N reaches 
100, and in the fourth case well before N reaches 150. 

Subsequent figures show the imaginary parts of the largest unstable 
eigenvalues of , for various 14(14) models (cf. (I5.0.2P ) as a function 
of 0. For all these truncations, we took N = 400. 

For comparison to our numerically constructed matrices, it is convenient 
to write out small block representations of our desired matrices. Given the 
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Plot of the rriax(lm(Eigeovalues)) vs. 


P3(V1)wittin = .2 


Plot ol the max(lm(Elgenvalues)) vs. N for P3(V2) 


. = 2.0 


oOgoooooooooooooooooooooooooooooo 


ooooooooooooooooooooooooooooooooott 


Figure 1: Left: The size of the imaginary part of an unstable eigenvalue, as 
a function of N, for an NxN matrix truncation for P 3 (Vi) with = .2. This 
is seen to stabilize strongly as soon as N is sufficiently large, say N > 100. 
Right: The size of such an imaginary part, as a function of N, for an 
NxN matrix truncation for ^ 3 ( 1 / 2 ) with 0 = 2. This is also seen to 
stabilize strongly as soon as N is sufficiently large, say N > 100. 
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Plot ol the (Im(Eigenvalues)) vs. N lor P4(V2) withn = 1.0 
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Figure 2: Left: The size of the imaginary parts of unstable eigenvalues as a 
function of N, for an x matrix truncation for ^ 4 (^ 1 ) with 11 = .1. This 
example has multiple non-real eigenvalues, but both are seen to stabilize 
strongly as soon as N is sufficiently large, say N > 100. Right: The size 
of the imaginary parts of unstable eigenvalues, as a function of N, for an 
NxN matrix truncation for ^ 4 ( 1 / 2 ) with H = 1. This is also seen to stabilize 
strongly as soon as N is sufficiently large, say N > 150. 
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^ 3 (^ 1 ) model, we have as the 6x6 truncation 


Mf = 


— A 3 


for 


bi = 


Cl = 


1 
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l{l + 2 ) 


(2^ + l)(2£-l) (2^ +l)(2£ + 3)’ 

and \i = l{i + \). Compare (I5.2.14p . There is a comparable construction 
for P‘i{V 2 ) using as in (I5.2.26P and (15.2.301) . 

For the F 4 (Vi) model, we take A{x^) = X‘i{x\ — 3/7), as in (I5.2.33p . As 
a result, for the operation of Mi on Q, we have 

^(3^3)0 = - ^2:3) Cr 

= (a£a£-i«r-2)C£-3 + + aj + a'j_^_i — -^^£-1 

+ fl£+l + o|+l + o|+2 “ -^^0+1 + («£+lO£+2a£+3)C£+3i 

(5.3.1) 

with ai as in p5.2.9p . and 


B{x 3)A Ci =- — {n +X4 A{x3))Ci. 

M 


(5.3.2) 
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Figure 3: Left: The size of the imaginary part for the the unstable eigen¬ 
value for P 3 (Fi), as a function of fl, with N = 400. Here, we see stability 
for n > .4, well below the Arnold stability bound. Right: The size of the 
imaginary part for the the unstable eigenvalue for P‘i(y 2 ) as a function of 
with N = 400. We see stability for > 2.25, which is again well below the 
Arnold stability bound. 


Again, for the reader’s convenience and to assist with comparison to 
numerically constructed matrices, we write down the 6x6 truncation for 
the P 4 (Vi) model: 


M? = 
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{2i+l){2l - 1) 

mi-2) 
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(f + l)(£-l) 


+ 


{i+2m 


{21 - l){2l - 2,) {2l + l){2l-l) {21 +‘i){2l + 1) 
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= 


1 i£ + l)i£-l) 

1 m-2) . 

1 

1 

1 {2£ + l){2£-l)y 

(22-l)(22-3)V 

1 

1 


and, again, + 1 ). 

For the -P 4 (i^) model, we also take ^(xs) as in (I5.2.33p . Then, we observe 
that for the operation of M 2 on Q, we have for multiplication by A[x'i) the 
same expression as in (15.3.ip but with 


an 


' (£ + 2 )(£- 2 ) 
(2£ + l)(2£-l)’ 


(5.3.3) 


as in (j5.2.26p and B(x 3 )A~^(£ as in ()5.3.2p . 

For the FliPVs) model, we also take A(x 3 ) as in (j5.2.33p . Then, we observe 
that for the operation of M 3 on (n, we have for multiplication by A(x 3 ) the 
same expression as in (j5.3.1jl but with 


an 


' (£ + 3)(£-3) 
(2£ + l)(2£-l)’ 


(5.3.4) 


and, again, B(x 3 )A~^(n as in (j5.3.2p . 

We turn to a discussion of the spectral results recorded in Figures [SHSl 
Figure El deals with the ^ 3 ( 14 ) models, for k = 1,2. In both cases, the graph 
indicates that M^ has a single eigenvalue with positive imaginary part, for 
12 in a certain range (0 < 12 < 0.4 for A: = 1, 0 < 12 < 2.25 for k = 2), 
and this imaginary part decreases monotonically to 0 as 12 runs over these 
intervals. In both cases, the imaginary part reaches 0 for 12 well below the 
threshold specified by the Arnold method, as worked out in ^4.11 Figure 0] 
deals with 24(14) models, for A: = 1, 2. In this case, we have two eigenvalues 
of Mfe with positive imaginary part, for 0 < 12 < 7.5 and 0 < 12 < 0.3 for 
A: = 1 and 0 < 12 < 3.4 and 0 < 12 < 7.0 for k = 2. Again, all imaginary 
parts reach 0 for 12 well below the threshold to which the Arnold stability 
analysis applies. In these figures, we see a new phenomenon. Namely, the 
imaginary parts of the eigenvalues are no longer monotone functions of 12. In 
fact, as seen in the bottom parts of Figured] (and also the top part of Fig. [5]) 
there is considerable oscillation of these imaginary parts, as a function of 
12, over certain ranges of 12, particularly for k = 2. The bottom part of 
Figure El has analogous graphs, for the 24(14) model, also illustrating such 
oscillation. 

It is our expectation that similar results hold for the operators M^, and 
hence F, arising in the linearizaton procedure of ^4.21 Going further, we 
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Plot of the (Im(Eigenvalues)) vs.n for P4(V1) witfi N = 400 



Plot of the (lm( Eigenvalues)) vs. Jlfor P4(V2) with N = 400 





Figure 4: Top Left: The size of the imaginary parts of the two largest 
unstable eigenvalues for P^iVi) as a function of with N = 400. Here, we 
see stability for H > 7.5, well below the Arnold stability bound. Top Right: 
The size of the imaginary parts of the two largest unstable eigenvalues for 
Piiy 2 ) as a function of fl. Here, we see stability for > 7, well below 
the Arnold stability bound. Bottom Left: A blow-up of the tail for the 
imaginary part of the largest unstable eigenvalue ^ 4 ( 1 / 2 ) to show that the 
oscillations towards stable are smooth at fine scales and not numerical errors. 
Bottom Right: A blow-up of the tail for the imaginary part of the smallest 
unstable eigenvalue Pi{V 2 ) to show that the oscillations towards stable are 
smooth at fine scales and not numerical errors. 
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Plot of tli8 (Im(Eigenvalues)) vs.n for P4(V3) ■ 


(ImfEigenvalues)) vs.il for P4(V3)' 



Figure 5: Left: The size of the imaginary part of the largest unstable 
eigenvalue for ^ 4 (^ 3 ) as a function of Here, we see stability for > 6.5, 
well below the Arnold stability bound. Right: The size of the imaginary 
part of the largest unstable eigenvalues for Pa^Vs) as a function of fl blown 
up to see that it is smooth curve. 


imagine there are further stability and instability results to be established 
for the Euler equation (jl.fl.ljl . We look forward to future progress on these 
problems. 
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Largest Imaginary Part of The Eigenvalues 


Plot of the max(lm(Eigenvalues)) vs. N for M1 (P3) withH = .2 





Largest Imaginary Parts of The Eigenvalues 


Plot of the (Im(Eigenvalues)) vs. N for M1(P4) withH = 0.1 







3St Imaginary Parts of The Eigenvalues 





Largest Imaginary Part of The Eigenvalues 


Plot of the max(lm(Eigenvalues)) vs. N for M2(P3) withH = 2.0 





Largest Imaginary Parts of The Eigenvalues 


Plot of the (Im(Eigenvalues)) vs. N for M2(P4) withH = 1.0 
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Plot of the (Im(Eigenvalues)) vs. Q. for M2(P4) with N = 250 
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